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INTRODUCTION 


The flight testing of airplanes in the past frequently 
has lent itself to methods of procedure somewhat less 
scientific than adventurous. It has always been a job 
fraught with uncertainties: the uncertainties of the 
weather, uncertainties of instruments and their func- 
tioning in flight conditions, uncertainties of scientific 
observation under conditions of discomfort and dis- 
traction, and uncertainties of the full meaning of the 
measurements made. The records obtained in flight 
have often been so inaccurate that no reliance could be 
placed upon them. The airplane may not have been 
flown steadily or long enough under steady conditions 
to establish a stabilized state of motion. The power of 
the engine has been indeterminate depending upon a 
great many variables, few of which, as a rule, have been 
measurable. The atmosphere is a fluid of uncertain 
motion and density. Drag items on the airplane, vari- 
ables in themselves, are frequently not under control. 
Finally, the technique of flight testing is still not very 
fully developed. 

Attempts made to compare the performance of two or 
more airplanes frequently lead to the conclusion that 
when the airplanes are not tested on the same type of 
day, with similar temperature conditions, the results 
are not comparable. Full-throttle ceiling may vary as 
much as fifteen per cent from day to day, and methods 
of reducing it to ‘‘standard conditions” are contro- 
versial. Level speeds have been so inaccurately meas- 
ured’and reduced to standard that the results of such 
testing at one time became usable only as very question- 
able sales figures. As a matter of fact, the airplane 
manufacturer himself frequently had no more than a 
vague idea of the real basic performance of his aircraft. 


Presented at the I.Ae.S. session of the A.A.A.S. Meeting, 
University of Washington, Seattle, June 21, 1940. 


HISTORICAL SKETCH OF FLIGHT TESTING 


One of the first attempts to reduce this flight test un- 
certainty was made by the British National Physical 
Laboratory. Early in the last war Tizard and Bair- 
stow developed the method used at the Experimental 
Aircraft Testing Station in England. There pilots 
were trained in technique and engineers were trained in 
methods of handling the flight data. At the same time 


flight test procedures were developed independently in 


France, and in 1918 there was a hot controversy over 
the superiority of the English vs. the French methods. 

In the United States two schools of thought de- 
veloped. The Army airplane testing activities early 
concentrated at Dayton, Ohio, where a group of test 
pilots and test engineers developed the Army method of 
flight testing and manipulation of the flight test data to 
give results obtainable under ideal standard conditions. 
This technique has been called the Army method of 
reduction. The Navy, starting from a somewhat theo- 
retical standpoint, developed both flight testing tech- 
nique and reduction methods on a different basis. 
This method has variously been called the Navy method 
or the Diehl method.* Efforts made to coordinate the 
Army and Navy methods were never very successful. 

All attempts to obtain greater precision in measuring 
aircraft performance gain their impetus largely from the 
still disconcerting discrepancy between flight test re- 
sults obtained by different pilots flying in different 
atmospheric conditions and reducing their data in 
different ways. Explanations for these differences and 
theories, to bring them into conformity, have been 
elaborate. Frequently the method used was to take a 
quantity of empirical flight data and from the per- 
centage variation determine constants which would 
reduce the discrepancy. 

Recent developments of the theoretical knowledge of 


* Engineering Aerodynamics by Walter Stewart Diehl; The 
Ronald Press Company, New York, 1936. 
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airplane performance variation with the factors influ- 
encing it have now provided a tool which, together with 
improved technique in obtaining accurate flight data, 
makes the measurement of airplane performance as 
accurate as other similar engineering measurements. 
The fact that it is not possible to climb as high with a 
given airplane on a hot summer day as on a cold winter 
day does not alter the basic performance of the airplane. 
PER- 


THE GENERALIZED CONCEPTION OF AIRPLANE 


FORMANCE TESTING 


Airplane velocity and climb-rate are the results of 
work done by the engine through the medium of the 
propeller. Since the level speed of the airplane is 
merely the condition of zero climb for any given operat- 
ing condition defined by the weight, power, propeller 
r.p.m., and density altitude, the performance of the air- 
plane for this operating condition can be expressed by a 
curve of climb-rate vs. velocity. Fig. 1 illustrates such 
a curve of climb-rate vs. velocity (C;, vs. V) for constant 
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Fic. 1. Climb and high speed performance. 


values of weight (W), power (P), propeller r.p.m. (JV), 
and atmospheric density ratio (c). A series of C, vs. V 
curves may be constructed to express completely the 
performance of the airplane for all possible operating 
conditions, each of which would be defined by a different 
set of the variables W, P, N, and co. 
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at sea level. 


The effect of weight and power on the performance of 
a typical airplane is illustrated in Fig. 2. The per- 
formance was computed for three combinations of 
weight and power for constant sea level density alti- 
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Fic. 3. Effect of altitude on performance. 
tudes and constant propeller r.p.m. for each power con- 
dition (but differing from one power to another). 

The effect of altitude (at standard conditions) is 
shown in Fig. 3 for the constant values of weight and 
propeller r.p.m., as given with the figure. 
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The points of zero = rate for the various operating 
conditions, as in Fig. 2, represent the respective veloci- 
ties in level flight. These points may be plotted on a 
chart of power vs. velocity for constant W, o, and N. 
Expressed in these terms the curve represents the 
fundamental of airplane performance. The curve of 
Fig. 4 shows the relationship of power and velocity in 
level flight for constant weight, W, constant relative 
density, o, and constant propeller r.p.m., NV. In order 
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required for level flight. 








AIRPLAN 


to show the effects of the variables, weight and density, 
upon the power required to produce a certain velocity, 
Fig. 4 has been replotted for a typical airplane in Fig. 5. 
Two weights and two typical altitudes are shown. 

The effect of the variable NV is shown in Fig. 6. This 
effect can be seen to be relatively minor and is due pri- 
marily to small changes in propeller efficiency. For the 
full range of densities and weights, a great many of 
these charts would be necessary were it not for the fact 
that the effect of weight and density can be eliminated 
by changing the factors power, velocity, r.p.m., and 
climb-rate, to parameters involving weight and density. 
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The first two of these parameters are P,,, and Vj, 
respectively, and are defined as follows: 
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where W = weight for which additional performance 
may be desired. 
W, =A standard reference gross weight for 


which the power-required curve is deter- 
mined for sea level density altitude. 

P = Brake power available. 

V = True velocity (true air speed). 


The fundamental performance curve, as given in 
Fig. 4, depends upon the three variables, W, o, and N. 
Changing the coordinates to P;,, and V;,, in accord- 
ance with the above equations, eliminates the variables 
W and o, and thereby effects a general applicability of 
the curve to all altitudes and all weights. The variable 
N must also be written in terms of W and oc. This 
introduces a third parameter, N;,, defined by the 


\ l/s 
relation N/a) . Fig. 7 may therefore be 


drawn, as in Fig. 4, to present a generalized power- 
required curve for a constant value of N;,,.. Variations 
of V;~ correspond exactly to variations in N of Fig. 6. 

The generalization effected by the introduction of the 
basic parameters involving density and weight can be 
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used also with the climb-rate relationship of Fig. 3. 
The climb-rate parameter involving these variables is 
termed C;,, and it is similarly formed by multiplying 
climb-rate by the square root of the density ratio and 
dividing by the square root of the weight ratio; that is, 


W\"'s 
Cry = 0 C;/ cz ) . Curves of climb-rate plotted 


in this manner are for constant values of P;, and Njy. 
Fig. 3 is replotted in Fig. 8 for constant values of Pj, 
and N;,. If the range of these curves is great enough, 
the complete plot will represent the total performance 
of the airplane. * 
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A further approach to fundamentals in determining 
airplane performance is possible by reducing weight to 
lift coefficient and reducing power required to drag 
coefficient. This is illustrated in Fig. 9. This familiar 
form of the performance curve is frequently referred to 
as the airplane’s “‘polar.”’ 

While this curve is not a true parabola, the part of the 
curve covering the operating range of the airplane so 


* The performance parameters appear simply and logically in 
the fundamental airplane performance equation for steady 
straight flight. (See N.A.C.A. Technical Report No. 408, 
General Formulas and Charts for the Calculation of Airplane Per- 
formance by W. Bailey Oswald, 1932. See also Applied Aero- 
dynamics by Leonard Bairstow; Longmans, Green and Co., 
New York, 1939 
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closely approximates a parabola that, for simplicity, the 
most useful part of the curve can be plotted as a 
straight line with the coordinates C;,* and Cp, as in 
Fig. 10. 

The “polar” of an airplane is not constant for the 
various conditions of power. Increase in power in 
general increases the drag of the portion of the airplane 
in the slipstream because of the increased air velocity. 
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The lift distribution is affected by the slipstream passing 
over portions of the wing. Change in the lift distribu- 
tion in general changes the induced drag. The pro- 
peller thrust in general has a lift component. The lift 
distribution and the lift component of the thrust affect 
the angle of attack necessary for a given lift coefficient. 
The turbulence and pulsating nature of the slipstream 
affect both the lift and drag. All of these power factors 
affect the polar. However, the first, the increase in 
drag directly due to increased air velocity in the slip- 
stream, is normally taken into account in the propeller 
efficiency. Such method is an approximation at the 
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best, but it is a simple one. If the profile drag coeffi- 
cient of the portion of the airplane included in the slip- 
stream is essentially constant, the slipstream velocity 
effect may be expressed as a constant factor of propeller 
efficiency. Cowl flaps and other variables complicate 
or invalidate this method. However, the use of this 
method defines propulsive efficiency, which is propeller 
efficiency reduced to account for the increase in drag 
directly due to increase in velocity in the slipstream. 
It may be defined as the overall efficiency of the pro- 
peller in developing thrust with respect to the entire 
airplane. Precisely this is not true because of the 
other effects ot the slipstream on the polar. When 
analyzing flight test data, these other effects are more 
simply represented as changes in the polar than as 
changes in propeller efficiency. Thus, a family of 
polars can be drawn for various power conditions, as in 
Fig. 11. 

Each polar of Fig. 11 is assumed to be for constant 
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power parameter (P;,,). They could equally well be 
drawn for constant values of the climb-rate parameter 
(Cj~). Variations of the r.p.m. parameter JN;,, intro- 
duce a family of curves for each of these conditions. 


THE CONCEPTION OF STANDARD ALTITUDE 


Standard altitude is defined by a sea level tempera- 
ture and pressure and a vertical temperature gradient. 
These are empirical figures supposed to represent a 
mean for temperate zone latitudes. These three quan- 
tities determine a definite relation between tempera- 
ture, pressure, density, and altitude. This relationship 
is illustrated in Fig. 12. For every pressure there is 
defined an altitude, and a scale of these altitudes can be 
used to calibrate a barometer. Then in the standard 
atmosphere this calibrated barometer may be used to 
measure altitude. In atmospheric conditions that are 
not standard, such a barometer offers an approximate 
measure of altitude. Altitude determined on this basis 
is called pressure altitude. 

Corresponding to this pressure scale there is also a 
density scale of altitude. If a density-measuring in- 
strument were available, it could be calibrated to 
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measure altitude accurately in a standard atmosphere. 
As in the case of the barometer, this might be used even 
in a non-standard atmosphere as an approximation. 
Altitudes measured on this scale are density altitudes. 

Every combination of pressure and temperature in 
the atmosphere defines a density and therefore a den- 
sity altitude. This suggests the plotting of a chart from 
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Fic. 12. Standard altitude characteristics. 


which density altitude can be read directly for any com- 
bination of pressure and temperature. Such a chart is 
illustrated in Fig. 13. 


THE REDUCTION OF FLIGHT TEST DATA 


A. Airspeed Calibration. The actual performance 
of the airplane may be measured in flight test, but this 
performance applies only to the conditions under which 
the test may be performed. The problem of reducing 
the flight test data is that of reducing the data to given 
conditions of weight, density, power, r.p.m., mixture, 
and any of the other several factors affecting perform- 
ance. 

It is necessary to correct for weight, because it is im- 
possible to perform any flight tests at constant weight, 
except performance in a truly power-off condition. 
The, operation of the engines continuously consumes 
fuel and thereby reduces the weight. The variation in 
temperature and pressure relation in the atmosphere 
from day to day, and even hour to hour or minute to 
minute, makes it impossible to complete any flight test 
under fixed temperature and pressure relations in the 
atmosphere. In other words, flight test performance 
data must be corrected for temperature and density. 


The power and r.p.m., of course, may be controlled 
within limits during the flight test; however, all of the 
possible conditions of power and r.p.m. cannot be cov- 
ered in any conceivably reasonable period of flight test 
time. This makes it necessary to interpolate flight test 
data and in some cases extrapolate the data for the de- 
sired conditions of power and r.p.m. Correction for 
engine operating mixture is also necessary. 

First of all, the instrument installation for measuring 
the air speed must be calibrated and corrected for the 
variables affecting it. The first of these variables is the 
weight. The angle of attack of the airplane for given 
indicated air speed, power, r.p.m., lift,* and drag* 
conditions is a function of weight. The greater the 
weight, the greater the angle of attack, and vice versa. 
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Fic. 13. Density altitude in terms of pressure altitude 
and temperature. 


The error introduced by the pitot installation normally 
may be considered purely a function of angle of attack. 
Of course it is conceivable that a pitot location might be 
such that other factors might also introduce appreciable 
errors. The usual location of the pitot in the vicinity 
of the leading edge of the wing toward the wing tip or 
the vicinity of the nose of the fuselage or hull is essen- 
tially free from these other errors. It is therefore clear 
that an accurate air speed installation calibration would 
be in the form of a plot of percentage error against angle 
of attack. Neglecting certain factors, the angle of at- 
tack is a function of the lift coefficient. In level flight 
the lift coefficient is equal to the weight of the airplane 
divided by the dynamic pressure and the wing area, 





* Lift and drag conditions may be varied by wing flaps, etc. 
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factor of one-half the standard sea level density is con- 
stant and for any given airplane the wing area is con- 
stant, at least for present normal designs. The lift 
coefficient then is reduced to a constant times the 


weight over the square of the true indicated air speed, 





{ Of the dynamic pressure the 
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Cr = |<) y3).* The calibration of the pitot in- 
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stallation may then be plotted in the form of per cent 
error against the factor W/V;*. An equivalent plot 
would be one of the ratio of the indicated air speed to 
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speed, —- This is illustrated in Fig. 14. Since the 
V; 

angle of attack for a given lift coefficient is affected by 

the deflection of the wing flaps, a separate calibration 


curve in general will be obtained for each flap deflection. 
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If the power has an appreciable effect on the lift for a 
given angle of attack, it will be necessary to develop air 
speed calibration curves for various conditions of power, 
including power-off. The introduction of the wing flap 
has considerably complicated the problem of air speed 
calibration. Besides that, the increasing appreciation 
for the effect of power has brought a realization of the 
added complication of power. The calibration, as 
plotted in Fig. 14, of course, is not directly usable. 
The usable form is the familiar one of the plot of indi- 
cated air speed vs. true indicated air speed, or, indicated 
+ For climb (or descent) a component of the weight is carried 
by the thrust (or drag), in which case the problem becomes more 
involved. For simplicity, the level case is here presented. The 
case of climb (or descent) is a simple extension of the case of level 
flight. 

* For the purposes of this paper, the true indicated air speed is 
defined on the basis of dynamic rather than impact pressure. 
This definition clearly has certain advantages. 
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correction vs. indicated air speed of pitot installation, 
and a complete calibration of this type will include a 
series of curves covering the range of weight, flap, and 
power conditions. Fig. 15 illustrates for a typical air- 
plane, a series of correction curves for power-off, for 
the extreme flap settings, and for various weights. 
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Fic. 15. Airspeed calibration chart 


B. General Performance. The data recorded during 
a flight test for performance include the basic factors of 
time, altitude, power, engine r.p.m., air temperature, 
carburetor mixture, and fuel load. These data may 
be recorded at regular intervals of time or altitude, or 
may be recorded continuously on a graph whose basic 
variable is time. In any case, the rate of climb is 
measured by the slope of the altitude vs. time curve. 
In general, this does not give a true geometric rate of 
climb. This is the case because the altitude, as meas- 
ured by the altimeter, is not the true altitude. Recall- 
ing the definition of standard atmosphere and the means 
by which the altitude calibration is made on the al- 
timeter, it is realized that the altimeter indicates true 
altitude only in case the pressure and temperature con- 
ditions result in a static pressure drop, below existing 
surface pressure, necessary to produce standard pressure 
at the particular altitude. However, so far as estab- 
lishing a definite figure for performance is concerned, 
it is only necessary that the standard temperature exist 
for the atmospheric pressure at which the performance 
is determined. In other words, each individual per- 
formance figure may be corrected for temperature at 
the particular pressure altitude. The airplane, pro- 
peller, and the engine are not at all interested in the 
geometric altitude. They are interested only in the 
temperature and pressure so far as the atmosphere is 
concerned. Neglecting for the moment the effect of 
moisture, air currents, etc., it can be shown that the true 
rate of climb is obtained from the rate of change of pres- 
sure altitude by multiplying by the ratio of the absolute 
temperature, existing during the test, to the standard 
absolute temperature, corresponding to the particular 
pressure altitude. This relation is easily obtained from 
the general gas equation and the hydro-static equation 
for the atmosphere. The air speed calibration deter- 
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mines the true indicated air speed. The pressure alti- 
tude and air temperature determine the relative density. 
The fuel load determines the gross weight. The engine 
gear ratio and the engine r.p.m. determine the propeller 
r.p.m. Combining these factors with the measured 
power and rate of climb, the basic parameters Vj, 
Piw, Niw, and Cj, are determined. Having deter- 
mined these parameters for all of the conditions tested, 
a plot similar to Fig. 8 results. Because of changes in 
weight due to fuel consumed, changes in temperature 
due to time and location, and changes in power due to 
ram, there is no series of points that may be directly 
faired to form curves at constant P;,, and N;,. The 
final curves that may be faired through such data result 
from a series of fairings, cross-plotting, and refairing, 
eliminating points widely in error, and applying judg- 
ment, realizing the accuracy that may be expected in 
flight testing. In order to illustrate this procedure, the 
successive steps will be described. 

First of all, it is possible to approximate curves for 
constant P;,, by the simple method of sketching in by 
eye as the surveyor does in drawing his contour maps. 
This gives the first approximation. On the basis of 
these curves, a cross-plot may be obtained, as shown in 
Fig. 16, plotting C;,, vs. Pj» at constant V;,, neglecting 
for the first approximation the effects of Nj». This 
cross-plot will at once reveal any gross error in contour 
sketching. These errors may be eliminated by a refair- 
ing of the original curves and a recross-plot. When the 
curves of Fig. 8 and the cross-plot of Fig. 16 are reason- 
able and consistent, a correction for N;. may be made. 
The details of this process become quite involved and 
their explanation would be too extended for inclusion in 
this paper. This correction for N;, roughly amounts 
to the following: The effect of variation of Nj» at 
constant P;, is a variation in propulsive efficiency. 
The variation of propulsive efficiency may be closely 
approximated by the use of propeller charts obtained 
from the wind tunnel tests on an installation very simi- 
lar to that tested in flight. In general, propulsive effi- 
ciency for the model tests may vary considerably from 
propulsive efficiency developed during the flight test. 
However, it is variation in propulsive efficiency with 
Ni» and not the actual magnitude of the propulsive 
efficiency that is desired. The variation in propulsive 
efficiency as obtained from the propeller charts may in 
general closely approximate the actual variation in 
flight. In any case, the effect of change of propulsive 
efficiency is small and any errors in calculating this 
effect, although even relatively large, say as much as 10 
or 20 per cent, will be negligible in the final result. 
However, in general the propulsive efficiency correction 
is in error much less than 20 per cent. 

Having corrected the points on Fig. 8 to certain con- 
stant values of N;,,, a refairing and cross-fairing process 
may be made to obtain the final faired results. This 
final fairing then represents the performance of the air- 
plane for all conditions of speed, power, r.p.m., altitude, 


temperature and weight within the range tested. 
Within reasonable limits, extrapolation may also be 
made. 
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Fic. 16. Cross-plot to fair performance curves of Fig. 8. 


C. Power Calibration. Once the data of Fig. 8 have 
been accurately determined, the performance for any 
given power may be easily determined. Whether any 
given power may actually be developed for given condi- 
tions depends upon the particular power plant. In 
order to determine whether the power plant will develop 
a certain power, a complete calibration in flight is neces- 
sary. The power rating of an engine, even though 
accurate for sea level and altitude conditions, does not 
in general apply to the installation in an airplane. For 
several reasons, the power output depends upon each 
particular installation. These reasons include car- 
buretor ram or its equivalent, engine cooling, cooling 
distribution, mixture (which may be different from test 
stand conditions in order to obtain desired engine cool- 
ing or economy), carburetor air intake temperature, 
exhaust backpressure, accessory loads, and any other 
minor factors. So far as actual power output is con- 
cerned, a power calibration is not necessary for part 
throttle. Part-throttle power calibrations will give 
power, r.p.m., manifold pressure, and mixture relations, 
but these relations are not necessary for the simple 
determination of performance which does not involve 
fuel consumption. Obviously performance relating to 
range, endurance, and economy necessarily must in- 
volve power plant fuel consumption calibrations. The 
basic calibration necessary for the determination of the 
performance limit of an airplane is the full-throttle 
power calibration for various conditions of r.p.m., alti- 
tude, mixture, indicated air speed, carburetor air tem- 
perature, and engine cooling. (Manifold pressure is a 
dependent variable which results from the calibration 
but is not necessary for its definition.) 

First, consider the power calibration at full throttle 
at best power mixture. In order to obtain this calibra- 
tion, it is necessary to obtain power at full throttle 
throughout the operating range of r.p.m. and altitude 
at the best power mixture. It is also necessary to ob- 
tain carburetor intake pressure as affected by ram for 
all of these conditions. The carburetor intake ram is 
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normally a function of true indicated air speed, weight, 
power, and density. The true indicated air speed is a 
major factor determining the pressure field about the 
airplane, and hence the pressure at the carburetor in- 
take. The true indicated air speed and the weight 
determine the angle of attack. The power output and 
the density determine the rate of intake airflow, and 
hence the pressure drop through the intake system. 
For a simple short-ducted carburetor intake, the effect 
of power on pressure drop may be negligible. The 
effect of weight depends upon the wing loading and the 
weight range. For the majority of designs this weight 
effect also may be essentially negligible. In case the 
power and weight effects are negligible, a single car- 
buretor intake ram curve may be developed as a func- 
tion of true indicated air speed. The intake ram cali- 
bration curve may be obtained and plotted as in Fig. 17. 
Another factor which may be applicable is that of flap 
deflection. The effect of flaps is to change the angle of 
attack for a given indicated air speed. Normally the 
condition of flaps deflected is not considered of impor- 
tance because prolonged operation with flaps deflected 


is unusual. 
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At full throttle the power output at a given pressure 
altitude is affected by the air temperature in two ways: 
First, the air temperature affects the density of the air 
entering the supercharger. The pressure rise through 
the supercharger is a function of this density. Hence, 
the pressure of the charge entering the engine is a func- 
tion of the density of the air entering the supercharger. 
Second, the density of the air entering the engine is 
also a function of the temperature. The power at full 
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throttle obtained under one condition of temperature at 
a given pressure altitude may be corrected to any other 
temperature by means of two corrections. The first 
involves a correction of the pressure rise through the 
supercharger and the second is a correction of the den- 
sity of the ingoing charge. Curves of the form of Fig 
18 have been developed for temperature correction of 
manifold pressure. Having obtained during the power 
calibration data on the relation between the density of 
the entering charge and the pressure rise through the 
supercharger, the pressure rise through the supercharger 
for any other intake temperature condition may be 
determined. It is therefore possible to obtain a series 
of curves at various engine r.p.m.’s, showing pressure 
rise through the supercharger as a function of carburetor 
inlet pressure at a given atmospheric temperature, as in 
Fig. 19. This serves as a pressure calibration of the 
supercharger. The power output of the engine at a 
given pressure of the charge entering the engine depends 
upon the temperature of that charge. This tempera- 
ture in general is a direct function of the temperature of 
the atmosphere. A relation which has been established 
for some time and which proves to be sufficiently accu- 
rate is that the power output of the engine is inversely 
proportional to the square root of the absolute tempera- 
ture of the atmosphere. 

The two temperature corrections, one for pressure 
rise through the supercharger and the other for the 
density of the charge entering the engine form the tem- 
perature power correction necessary for the power cali- 
bration. 

Because the air intake system of a power plant is of 
necessity close to a source of heat, the temperature of 
the air may be raised above that of the free air. Any 
rise in air temperature from this cause will cause loss in 
power as compared to that which would be developed 
were there no rise. However, such temperature rise 
normally is small and does not affect the method of 
power calibration. It appears only as a defect in the 
power plant installation that acts as one of several fac- 
tors to reduce the power output as compared to test 
stand calibration. * 

Because of the latent heat absorbed by the rapidly 
vaporizing fuel, the region after the carburetor is an 
excellent refrigerator under certain circumstances of 
carburetor design, power, throttle, and outside air 
temperature. If the air carries sufficient moisture as 
vapor or moisture in proper quantity as droplets, ice 
may be deposited on the duct walls. The ice deposit 
may grow to such an extent that it completely chokes off 


* This method differs from that in Aircraft Engine Power 
Estimation from Intake Manifold Density by Rhines, Thomas B., 
and Schipper, Peter W., Journal of the Aeronautical Sciences, 
Vol. 6, No. 3, page 106, January, 1939, because of its purpose. 
Its purpose is to develop an engine calibration easily usable with 
airplane flight test performance reduction. The purpose of the 
method of Rhines and Schipper is to develop a basic method of 
study of engine power at altitude usable basically for engine per- 


formance. 
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Fic. 18. Manifold pressure correction for carburetor air 
temperature. 


the intake. This ice formation must be prevented in 
some manner. The usual method is to provide a means 
of heating the intake air. Generally such a provision 
not only raises the temperature but acts as a choke and 
reduces the intake pressure. For some carburetor in- 
stallations for which preheat operation is not unusual, a 
power calibration with the preheat on may be necessary. 
Basically this is no different from the normal calibration 
for cold air. 

Another power correction which is necessary in ob- 
taining a complete power calibration is that for mixture 
of the intake charge. Normally the primary power 
calibration of an engine is made at the mixture for best 
power. For mixtures appreciably richer or leaner than 
best power mixture, the power is measurably reduced. 
The effect of mixture on power is essentially independ- 
ent of altitude. It is also independent of temperature 
so long as the detonation limit is not reached. Hence, 
normally it is sufficient to obtain the effect of mixture 
on power throughout the operating range of power and 
r.p.m. without particular regard to temperature or alti- 
tude. The plot of Fig. 20 illustrates the effect of fuel/ 
air mixture on power.{ If the detonation limit is 
reached the leanest mixture is determined by the atmos- 
pheric temperature. 

Having established the principles by which the power 

t Good accuracy in measuring fuel/air ratio is obtained if the 
instrument is serviced and calibrated before and after each flight. 
Otherwise, fuel/air ratio measurements may be of very doubtful 
accuracy, 
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Fic. 19. Pressure calibration of engine supercharger. 


calibration is obtained, it would be of value to follow 
the steps of obtaining actual power calibrations. First, 
of course, the power data must be obtained in flight 
tests for the atmospheric conditions existing at the 
particular time. From this data, power curves may be 
developed for any of the series of temperature condi- 
tions in the atmosphere. The actual power obtained at 
any particular point may be corrected to a desired indi- 
cated air speed by means of the intake ram calibration 
and the supercharging pressure rise calibration. The 
use of these two calibrations gives the correction to the 
pressure of the charge entering the engine. How the 
power is affected by the pressure of the charge is deter- 
mined from the relation between so-called manifold 
pressure and power. This relation may be determined 
either by tests on the ground or by tests in the air 
Normally, the first is available from the engine manu- 
facturer’s power calibration of the engine and the second 
is available from the power calibration in flight. One 
offers a check against the other. The relation between 
manifold pressure and power is illustrated in Fig. 21. 
Now, having determined the power output at a given 
indicated air speed, the next step is to obtain the power 
at a given air temperature. First, a correction for 
manifold pressure is necessary and is obtained by means 
of the curves of Fig. 18. Second, a correction for den- 
sity of the charge entering the engine is necessary and is 
obtained by means of the relation that power is inversely 
proportional to the square root of the absolute tempera- 
ture of the atmosphere. Having corrected a series of 
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points in this way at various altitudes at full-throttle 
and at a given r.p.m., it is possible to plot a curve of full- 
throttle power vs. altitude for a given condition of indi- 
cated air speed, r.p.m., and atmospheric temperature. 
A series of these curves may be obtained for the range of 
operation. Such a plot is illustrated in Fig. 22. Hav- 
ing obtained such a complete power calibration, it is a 
simple process to determine the performance of an air- 
plane by means of the charts of Fig. 8 and Fig. 22. 
Thus, by means of two such complete charts, the entire 
performance of an airplane may be presented for opera- 
tion at best power mixture. For operation at mixtures 
other than best power, a correction for power by means 
of a series of curves, as in Fig. 20, will be necessary. 
Two primary curves supplemented by the mixture 
curve then present the full range of performance of the 
airplane for normal flight. A number of factors have 
been neglected in reducing the data in this manner. 
Among these factors are Reynolds Number, the effects 
of compressibility on the aerodynamics of both the air- 
plane and the propeller, the variation in power output 
from engine to engine, the variation in power output 
during the life of the engine, and other such factors. 
The use of the torquemeter will eliminate the power 
variation problem. 


THE TECHNIQUE OF HANDLING THE AIRPLANE FOR PER- 
FORMANCE TESTS 


Few other fields of engineering measurement offer so 
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many possibilities of error as the determination of air- 
plane speed and climb-rate. In the past the precision 
of such tests has sometimes been so bad that results 
could not be relied upon within 10 per cent because of 
lack of weight control, improper instrumentation, un- 
controlled technique, or failure to evaluate the variables 
involved in the measurement of performance. When 
modern methods and technique are available a flight 
test accuracy of +1 per cent can be realized in measur- 
ing both climb-rate and velocity. 

The principal factors requiring careful control in flight 
testing are (1) weight control such that the gross weight 
of the airplane can be determined at any time during the 
test, (2) instrumentation provided in advance so as to 
give accurate measurements of every element vitally 
affecting the space-energy picture, (3) drag control so 
that parasite resistance does not change in some un- 
known manner during the test, (4) restriction of testing 
to known atmospheric conditions, particularly eliminat- 
ing vertical air currents, (5) maintenance of steady con- 
ditions of flight long enough to assure stabilized energy 
relationships, and (6) experience in the type of precision 
flying required for testing. 

Weight control requires accurate weighing and subse- 
quent checking of all items added or subtracted from 
the airplane. Primary on all large aircraft is the item 
of fuel used during the flight. The difference in weight 
between the start and end of a flight test will often be 
enough, if uncorrected, to make otherwise excellent 
results inconsistent and unreliable. On each flight, a 
data sheet is maintained for the purpose of keeping a 
running time-log of the gross weight of the airplane. 
Fuel consumed can be read from fuel gages from time to 
time, and weight computed from these figures and 
checked at the end of the flight. Such a log and its 
accompanying chart of weight vs. time are shown in 
Fig. 23. Each reading of velocity and climb can 
then be corrected for the proper ratio of gross weight to 
standard gross weight. 

Instrumentation for tests requiring accuracy can be 
done properly only when it is begun prior to the begin- 
ning of the test, preferably prior to completion of the 
airplane. The principal performance instruments re- 
quiring careful selection, calibration, installation, and 
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attention are (1) the air temperature gage, (2) the sensi- 
tive altimeter, (3) the air speed indicator, (4) the tachome- 
ter, (5) the torquemeter, (orin its absence the mani- 
fold pressure gage, the carburetor air temperature gage, 
the mixture indicator and engine temperature indica- 
tors). 

Atmospheric temperature is a very difficult measure- 
ment to make accurately. Small errors in this meas- 
urement are equivalent to excessively large errors in 
density altitude and air speed. In order to show this 
effect, the following table has been made: 

Error Introduced by 
Common Error Temperature Error 


Air Temp.—=10°C. for improper 


measurement 

Air Press.— = 20 ft. on std. sensitive + 1200 ft. 
alt. 

Air Speed +1 m.p.h. on test instru- +3.5m.p.h. at 200 m.p.h. 
ment 


Common errors in obtaining air temperature accu- 
rately introduce errors in altitude 60 times the com- 
monly obtained accuracy of making this reading and 
31/2 times the commonly obtained accuracy in obtaining 
air speed. 

It is clear from this that the principal effort should be 
placed upon reducing the usual error of temperature 
measurement by a large amount in order to equal the 
accuracy of other important measurements. An ac- 
curacy of +1°C. is desired at the present state of the 
art. Indicated temperature differences as great as 
35°C. have been noted on an airplane in flight on differ- 
ent thermometers whose elements were located at vari- 
ous points on the airplane. 

There are several primary reasons for the excessive 
error in obtaining air temperature. One of these is the 
location of the thermometer element at a point where it 
will not measure the true air temperature at the altitude 
of the test. The element may be picking up heat from 
the engine or from the slipstream of the propeller or 
from the sunlight (either direct radiation or conduction 
of heat from surrounding sun-heated surfaces) or from 
adiabatic temperature changes owing to compression or 
expansion at localized points around the airplane. The 
avoidance of the engine, propeller, and sunlight effects is 
readily obvious. Elimination of the errors due to 
adiabatic heating can be accomplished in one of several 
ways. Since compression effects upon temperature are 
virtually non-existent at speeds under 100 m.p.h., the 
airplane can be flown slowly at the altitude in question 
long enough to get stabilized temperature readings. 
An experimental correction can easily be determined 
which can be used for this purpose if the temperature 
element is normal to the airstream. Some test per- 
sonnel prefer to make the temperature reading from a 
separate slow flying airplane sent up for the purpose of 
obtaining the altitude-temperature gradient at the time 
of the test. 


Lag effects are quite serious in climbs and descents 
where the temperature indications may be several de- 
grees behind the actual temperature at any instant dur- 
ing the changes in altitude. These effects can some- 
times be satisfactorily reduced on slowly climbing air- 
planes by measuring temperatures during both climbs 
and descents made at the same rate of change of alti- 
tude. 

Thermocouples have been found to be more satis- 
factory than resistance type elements for measuring 
outside air temperature owing to the difficulty of cor- 
recting the latter for the effects of changes in rates of 
radiation in air of varying density. When this element 
is normal to the airstream it is subjected to total pres- 
sure developed by the velocity of the airplane. 

When so installed, a simple experimental correction 
can be made for the rise in temperature due to compres- 
sion. If the element is not normal to the airstream, the 
correction becomes involved by local pressures. 

Accurate pressure instruments for measuring pressure 
altitude are now available so that this measurement 
should be the easiest of all to obtain within the desired 
accuracy. The standard cockpit sensitive altimeter is 
perfectly satisfactory for determining pressures if prop- 
erly calibrated and installed. It must, of course, be 
vented to the static tube of the air-speed pitot in order 
to establish a standard. Of course this standard must 
be calibrated in the same manner as the air speed to 
eliminate errors due to local pressures. 

A special air-speed indicator is desirable for perform- 
ance testing primarily for convenience and saving of 
flight test time. Graduations in one m.p.h. intervals 
are desired. The calibration of this instrument is made 
separately in addition to the calibration of the complete 
air speed indicator installation. Then, when a new 
instrument is used for replacement, it will be unneces- 
sary to recalibrate the airplane. 

The simplest form of air-speed indicator calibration is 
made by flying the airplane at its full gross weight, at 
very low altitude, across a measured course at constant 
indicated air speed, and timing the runs for a distance 
great enough to reduce the error of measurement to a 
value consistent with the error of other flight test meas- 
urements. When the airplane is flown across such a 
course both up-wind and down-wind a sufficient number 
of times throughout the speed range required for the 
calibration, a curve can be drawn of indicated air speed 
vs. true indicated air speed. With the usual instru- 
ments, such a method of calibrating the air speed indi- 
cator involves, however, so many and such large errors 
as a rule that it is not feasible for accurate flight testing 
unless ideal conditions are available. An automatic 
timer of high accuracy allows the use of this method 
over a short course of a few hundred feet. In such a 
short distance variation in some of the factors is mini- 
mized. But even so, the wind conditions must be very 
consistent or constant. 

One modern method of air-speed meter calibration is 
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the use of a trailing bomb. A trailing bomb is a pitot 
tube mounted in such a manner in a streamlined, finned 
body, that it can be trailed below and behind the air- 
plane in a region undisturbed by the airflow around the 
wings or body. The disturbance has been found to be 
less severe behind and below than directly below the 
airplane. This fact fits well with the trailing position 
of the bomb hung by a cable from the rear of the air- 
plane body. The cable carries the tubes transmitting 
the pressures from the pitot to the airplane air speed 
indicator. Cloth streamers are tied at intervals to the 
cable to give the desired trail angle and bomb position 
with respect to the airplane. A cable length at least 
equal to one span should be used, and the bomb should 
be lowered from a point as far aft as possible. 

The calibration of an air-speed meter installation 
with a trailing bomb may be accomplished at any alti- 
tude. The bomb is first calibrated in a wind tunnel, 
from which a relation between indicated air speed and 
true air speed is determined. 

The calibration technique consists in obtaining at a 
series of air speeds throughout the speed range simul- 
taneous readings of two air speed indicators, one of 
which is connected to the airplane’s pitot-static and 
the other connected to the pitot-static of the trailing 
bomb pitot. Care is necessary to avoid trailing the 
bomb in air of velocity different from that at the air- 
plane level. Also, unless the airplane and bomb pitot 
and static lines are properly balanced for changes of 
pressure, the trailing bomb air-speed calibration must be 
performed in level flight. 

Of primary importance in measuring performance, 
both climb-rate and speed, are the instruments dealing 
with the power available at various altitudes and atti- 
tudes of the airplane. The tachometer is now suffi- 
ciently developed and standardized as to give results 
well within desired flight test accuracy. There are, 
however, a number of tachometers on the market whose 
accuracy should be accepted with caution. Until 
torquemeters came into use for flight testing, measure- 
ments of torque had to be made indirectly from readings 
of intake manifold pressure, carburetor air temperature, 
fuel mixture strength and other variables entering into 
the determination of b.m.e.p. The methods of ac- 
complishing this are explained above in the section 
dealing with engine calibration. 

The instruments themselves offer no serious difficulty 
from the viewpoint of calibration and interpretation. 
All instruments should be re-calibrated at intervals in 
order to assure their accuracy. An able man assigned 
to the instrument problem can easily save large losses 
arising from uncontrolled instrumentation. 

The problem of drag control should be so obvious as 
to require nocomment. It is, however, a problem to be 
watched carefully. The best means of- watching it 
during development tests of a new airplane has been 
found to document each flight very completely from 
the viewpoint of changes occurring to the external char- 
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acteristics of the airplane. During flight notations may 
be made frequently of cowl flap position, landing gear 
position, or any other variable which might affect flight 
characteristics. 

The atmosphere is a sea of air in which man has, after 
many years, managed to navigate unsupported by solid 
connection to the earth or by means of any appreciable 
buoyancy of simple displacement. This sea is not 
quiescent but is ever in a turmoil under the forces of 
thermo- and aero-dynamics that visibly produce our 
weather. It may flow as a tumbling waterfall over the 
mountainous areas or become a turbulent, rolling, seeth- 
ing mass, or flow smoothly and quietly with unappreci- 
able ripple. The air currents may be made visible by 
the clouds they produce, or by the smoke, dust, soil, and 
sand that they lift and carry. Even if visible, the mag- 
nitude of their effect is not directly apparent. Visible 
or not visible, the effect of air currents has been fre- 
quently ignored in determining airplane performance. 

The effect of air currents upon the accuracy of per- 
formance tests can be figuratively disastrous. The 
way in which air currents affect the performance of an 
airplane is well illustrated by the feats of soaring. We 
might well investigate some quantitative figures of air 
current effects. 

(1) Vertical Currents. Vertical currents are those by 
which thermal soaring is possible. A vertical current 
adds directly, vectorially, to the rate of climb. In 
level flight the speed is changed by that amount neces- 
sary to produce a vectorial sum of zero. For example, 
in Fig. 2, one may quickly determine the change in 
level speed for a descending current of 500 ft./min. by 
finding the difference between normal level speed and 
the speed for 500 ft./min. rate of climb for the given 
conditions of power, weight, and altitude. 

A vertical current may have a horizontal gradient. 
Its effect may be appreciable. 


W_ dw ; ; ; 
—-xX—-xXVXV=Wx AC 
ds 
AC_ Vy dw 
Vg’ ds 


where 


W = airplane weight 
g = acceleration due to gravity 


w = vertical current 
s = horizontal distance 
V = air speed 


A C = change in rate of climb 


If the horizontal gradient is 1.61 ft./sec. (about 100 ft., 
min.) per 1000 ft. for an airplane flying at 200 ft./sec. 
(about 140 m.p.h.), the rate of climb is increased by | 
per cent of the air speed or 120 ft./min. 

(2) Horizontal Current. <A horizontal current of con- 
stant magnitude with respect to all axes adds simply 
vectorially to the ground speed but has no other effect. 
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If a horizontal air current has a horizontal accelera- 
tion of given value, the airplane performance may be 
measurably affected. Taking the wind velocity and 
airplane velocity positive in the same direction, a posi- 
tive wind acceleration will decrease the airplane per- 
formance. This is true since the airplane performance 
depends on its momentum with respect to the surround- 
ing air. If the air current is continually increasing, the 
airplane finds itself struggling to maintain its momen- 
tum with respect to the air. In other words, it spends 
part of its thrust in accelerating with respect to the 
ground. This required acceleration might be expressed 
in equivalent rate of climb. 


h/ 
: XaXV=WXAC 


Ss 


AC=a‘*y 
g 


where 
a = acceleration of air current 


A horizontal acceleration of 1/10g changes the rate of 
climb of an airplane flying at 180 ft./sec. (about 120 
m.p.h.) by 1080 ft./min. A horizontal acceleration 
might easily give the impression of a strong vertical 
current. 

A vertical gradient of a horizontal current does not 
affect level air speed but may very appreciably affect 
rate of climb. An airplane climbing in a current of such 
gradient is in the same position of an airplane flying in 
an accelerating current of zero gradient. If the vertical 
gradient is positive, the airplane has to exert itself to 
maintain its air speed. 


Wd 
~xD2xcxV=WKaAc 
g dh 
ac lye 
Cg dh 
where 
h = vertical distance 
C = rate of climb in still air 


For a vertical gradient of 16.1 ft./sec. per 1000 ft. for an 
airplane flying at 200 ft./sec. (about 140 m.p.h.) the 
rate of climb is reduced by 10 per cent. 

An angle of yaw on the air-speed pitot may introduce 
quite an appreciable error into the air-speed indicator 
readings. When the air-speed pitot is located in such a 
position relative to the wing or body that flight at an 
angle.of yaw influences the pressure distribution in the 
region of the pitot head, the error introduced by yaw is 
even greater. Errors of +4 m.p.h. are not uncommon 
for barely appreciable angles of yaw on some aircraft. 
For steady level flight, the air speed for a stabilized 
zero-yaw condition will be different from that at a 
yawed-flight condition, primarily because of the fact 
that the airplane is sideslipping and must therefore be 


held in a climbing attitude in order to maintain the 
same altitude. For all of these reasons, it is quite 
necessary to maintain zero angle of yaw. This can be 
done either with a yaw meter, or for essentially sym- 
metrical conditions, by means of the standard turn and 
bank indicator or artificial horizon. 

Unless the test pilot has a definite knowledge of the 
precision necessary in taking readings of his various 
instruments, he cannot arrive at a satisfactory precision 
with the most economical use of his time. It is not 
always possible to wait until the atmospheric conditions 
are perfect before making a flight test. It is also not 
always possible to throw out flight tests in which verti- 
cal currents were encountered. The great cost of flight 
test time requires that the test pilot determine the 
minimum time required to obtain the desired accuracy. 
For level flight readings, it should be possible to obtain 
an accurate reading of + one-half per cent. It should 
be possible to repeat flight tests and get results with this 
degree of precision. Due, however, to the large number 
of variables which cannot always be brought under 
control from day to day, any comparison of speeds 
made for the purpose of obtaining drag indications must 
be made over a full range of power output. This is 
especially true when torquemeters are not available. 

Due to the very large differences in lag and sensitivity 
of the cockpit instruments, it is essential for accurate 
performance measurements to guide the technique of 
piloting by means of the proper sequences of instrument 
reading. It is frequently impossible, for instance, to 
fly in any given steady attitude by means of the rate-of- 
climb meter because of its lag and enormous variability 
with airplane angular velocity. It is virtually im- 
possible for a pilot to guide a transport airplanes’ de- 
scent by the rate-of-climb meter alone. It is necessary 
in this case to determine the approximate air speed 
which gives a given average rate of descent, and then 
change this air speed in accordance with the average 
indications of the rate-of-climb meter by a trial and 
error process until the desired rate of descent is ob- 
tained. 

This technique can be acquired only by practice in 
obtaining the desired stabilized conditions. It is for 
this reason that testing of performance with a pilot not 
trained in test flying is a very wasteful process from the 
point of view of flight test time and flight test accuracy. 


THE TESTING OF CLIMB PERFORMANCE 


The flight tests for determining climb performance 
require the control of weight throughout the flight, as 
explained above, so that the instantaneous weight can 
be used for climb-rate reduction. A continuous record 
of fuel load is necessary in order to prepare the plot of 
gross weight vs. time from calculations which are made 
to obtain the coordinates C;,, and V;,,.. The layout of 
the flight plan is based upon the principle that optimum 
climb-rate may be determined by a series of short steady 
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Fic. 24. Observed climb-rate for various constant indi- 
cated airspeeds. 


climbs at predetermined air speeds and powers. The 
air speed chosen shall bracket the speed for best climb- 
rate. The chart of altitude vs. time for such a test will 
have the appearance of a sawtooth, as in Fig. 24. The 
slope of the climb portions of the record give the rates 
of climb at the altitude past which the steady climb 
attitude was maintained. Frequently it is possible to 
obtain greater accuracy than the barograph record can 
give by having an observer make a continuous plot of 
altitudes at 15-second intervals. The standard sensi- 
tive altimeter is more accurate than the usual flight test 
barograph. Its use also offers the advantage of ena- 
bling the observer to discard portions of the climb obvi- 
ously affected by vertical currents or by inaccurately 
maintained velocity observable while plotting results. 

Constant power is maintained in climbs made below 
critical altitude or at partial powers to be used in the 
density climb-rate curve of Fig. 16. Sawtooth climbs 
are required at several altitudes and at a sufficient num- 
ber of air speeds in order to insure (a) that best climb- 
rate speed has been bracketed and (b) that an accu- 
rately formed curve of velocity vs. climb-rate is obtained 
in the region of best speed. Sawtooth climbs at re- 
duced power are necessary in order to determine the 
full range of climbing ability of the airplane. 

On multi-engined aircraft, climbing performance for 
one or two engines inoperative is determined by making 
sawtooth climbs under the desired conditions. For an 
accurate determination of ceiling with one or more 
engines inoperative, sawtooth descents above ceiling 
are desirable, and rates of descent under the steady 
conditions are determined in the same manner as the 
rates of climb were determined. 

The so-called check climb, frequently made to check 
upon the accuracy of the sawtooth climb-rate deter- 
minations, is actually an unnecessary flight since it 
determines, as a rule, not the climbing ability of the air- 
plane, but the average accuracy of the pilot who makes 
it. Comparisons between airplanes made by comparing 
check climbs made by different pilots are seldom accu- 
rate engineering measurements. 

Below critical altitude the “rated-power’’ climb will 
have a similar shape to the constant-power climbs 
throughout the density-altitude climb-rate chart. 
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Cruising-power climb is a characteristic curve that is 
very useful for all airline work. In fact, rate-of-climb 
vs. density-altitude curves for various cruising powers 
are needed in order to complete the picture of airline 
performance. Full-throttle rate-of-climb curves can 
be drawn for various temperatures above and below 
standard. This will complete the picture of the full 
range of climb performance. 

When best angle of climb as well as best rate of climb 
is desired, as it is for take-off power, rates of climb 
should be determined at several velocities less than the 
velocity for best climb-rate. The curve of climb-rate 
vs. velocity can be extended downward to the point of 
tangency of a line drawn from the origin. This point of 
tangency defines the speed for best angle of climb. 
This is primarily important only for take-off power, 
where obstacle clearing ability is to be determined. 

It may also be important to have rates of climb deter- 
mined accurately for speeds higher than the speed of 
best climb-rate. These rates of climb are useful in 
determining airline operating flight paths where 
cruising climbs are carried out at higher speeds in order 
to get the maximum distance covered in a given length 
of time for the first part of the airline flight. 

The traditional method of determining rates of climb 
from the power-required curve and the power available 
relationship is illustrated in Fig. 25. Excess power here 
assumes a static condition of the power-required curve. 

As was pointed out in the section on the generalized 
conception of performance, this conception of ‘‘excess 
power” is an approximation. The power-required 
curve is not representative for all conditions of power; 
generally it is the level-flight power-required curve that 
is determined. Accurately, there will be a series of 
power-required curves throughout the range of applied 
excess powers. 
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An important factor in making the power required for 
climb differ from that required for level flight is that at 
constant indicated air speed in a climb, the true air 
speed is increasing. This increase in air speed is ac- 
celeration that requires thrust to produce. Therefore, 
in a climb at constant indicated air speed (or at any 
changing indicated air speed that corresponds to chang- 
ing true air speed), additional power is required for 
acceleration. 


VELOCITY IN LEVEL FLIGHT AT VARIOUS POWERS 


A determination of level flight speeds, at various alti- 
tudes, powers, gross weights, propeller r.p.m., engine 
b.m.e.p., and fuel flows, is a very important part of the 
flight test program. In it the principal correlation 
between wind tunnel and flight test is determined. If 
the airplane is equipped with torquemeters, this deter- 
mination is greatly simplified, and the accuracy is 
greatly improved. It requires, nevertheless, very 
careful instrumentation and technique. It also re- 
quires careful control of weight, since there is fre- 
quently a very large difference in the level speed at any 
altitude for different gross weights. Cruising speed is, 
of course, to be plotted against density altitude, and the 
simplest form of such a plot is shown in Fig. 26 where 
pressure altitude and atmospheric temperature deter- 
mine the density altitude for use throughout the velocity 
chart. On such a density altitude basis, velocities for 
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Fic. 26. Cruising power chart 45,000 lbs. gross weight. 


constant power can be accurately plotted up to the full- 
throttle conditions. Full throttle will not occur at the 
same point in this chart if the atmospheric temperature 
varies from standard. It is, however, possible to draw 
full-throttle lines for constant variations from standard 
temperature. 

On airplanes equipped with fixed-pitch propellers, or 
with two-position controllable propellers, or with pro- 
pellers having a pitch-angle indicator, cruising charts 
can be drawn showing the velocity relationships for con- 
stant propeller r.p.m. and constant power, and at the 
desired propeller blade angle. The propeller r.p.m. 
lines will appear, in this case, as in Fig. 27. This forms 
a very satisfactory method of cruising control on small 
airplanes, or on those without elaborate instruments. 
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Fic. 27. Cruising power chart for airplane equipped with 
fixed-pitch propellers 


For airplanes having constant-speed propellers, it is 
preferable to draw cruising charts with constant values 
of b.m.e.p. since the speed can be altered to give the de- 
sired b.m.e.p. for any power. The advantage of oper- 
ating at constant b.m.e.p. is that in this condition the 
operating severity on the engine remains at a fixed de- 
terminable value. Another advantage is that specific 
fuel consumption can be determined quite accurately 
and maintained at any desired value. On such a chart 
it is possible to draw lines of constant power which are 
also lines of constant r.p.m. and fuel flow. From these 
lines it is simple to determine range, and lines of con- 
stant miles per gallon can then be superimposed upon 
this chart. This is illustrated in Fig. 28. 

The use of charts of the nature of Figs. 26, 27, and 28 
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Fic. 28. Cruising fuel consumption and range chart. 


for the purpose of cruising control, requires torque- 
meters, power calculations, or air-speed correction cal- 
culations, or their combination, depending on the ac- 
curacy desired or justified. 

It may be well to note that in Fig. 26 there is a sharp 
change in the nature of the curves of constant power 
above the full throttle line. This is caused by the 
change in propulsive efficiency due to the increase in 
r.p.m. necessary to maintain constant power. 


DISCUSSION 


The situation, so far as flight-test performance, flight 
testing, and reduction are concerned, is not so well 
settled as might appear by the facts presented in this 
paper. The methods presented here have been de- 
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veloped only for the simple gear-driven supercharger. 
For some time development work has been progressing 
on exhaust-driven turbo superchargers. The factors 
that affect the power output of such an installation are 
more numerous and more involved than for the gear- 
driven supercharger. Some excellent testing has been 
done on the ground. Other testing has been done and 
is being done in the air. The evaluation of all of the fac- 
tors involved has not yet been determined and will be 
the basis of extensive research in the future. 

The use of the hydraulic drive for remote drive and 
for gear reduction is now on the horizon. Its regular 
service use appears not too farin the future. For pres- 
ent engines there is a direct relation between the engine 
speed and the propeller speed. The introduction of the 
hydraulic drive will eliminate this simple situation. 
Another factor of complication will be introduced into 
the reduction and interpretation of airplane perform- 
ance. 

The variable area wing is now being used. 
developments along this line will go is not immediately 
apparent. The use of any normal wing flap, and in 
particular the Fowler flap, allows essentially a change in 
the wing area. Therefore, the performance of an air- 
plane with a variable area wing will be handled in the 
same way as the performance for various wing flap set- 


How far 


tings. 

The constant trend toward increasing speeds has al- 
ready brought airplane performance into the outer 
fringes of the effects of compressibility. The entire de- 
velopment of the methods presented in this paper has 
essentially neglected the effect of compressibility of the 
atmosphere, except insofar as minor corrections may 
account therefor. In other words, the atmosphere has 
been treated essentially as an incompressible fluid. 
When the local velocity of any part of the airplane or 
propeller reaches the speed of sound for the air at the 
particular locality, the air flow takes on a sharp change 
in character, and it results in extremely great changes in 
lift and drag. This compressibility effect may be 
avoided by proper aerodynamic design up to a certain 
point. As aerodynamic design resources are exhausted 
and airplane speeds increase, it eventually becomes im- 
possible to completely avoid the effects of compressi- 


bility on all parts of the airplane. These effects of com- 
pressibility do not lend themselves to the simple laws 
that have been used in the development of the methods 
in this paper. The complete answer to the question on 
how this effect will be handled in the interpretation and 
reduction of airplane performance does not appear im- 
mediately. 

The method of obtaining, presenting, and reducing 
flight test performance data at altitude that has been 
covered by this paper gives a means of checking the 
accuracy of each point against all of the other points. 
The data obtained at any given altitude are not con- 
sidered as data separate and by themselves. They are 
considered in relation to all of the data obtained at all 
of the altitudes tested. By this means the inaccuracies 
of any individual point are reduced toa minimum. A 
principle that may be used in preparing the data, one 
that is the basis of the scientific method, is that of 
checking the data against theory. By checking the 
theory against the flight test and flight test against the 
theory, it is possible to obtain a rational shape to the 
curves that are faired to the data. Having obtained a 
reduction of the data that is consistent with theory, it is 
possible to interpolate very accurately for any condition 
of weight, power, altitude, and temperature and ex- 
trapolate with good accuracy and to within reasonable 
limits. Thus, the complete performance of the airplane 
at all weights, powers, altitude, and temperatures may 
be obtained by means of specific tests at particular alti- 
tudes covering the operating range. The variations in 
performance with weight, power, and temperature are 
obtained on the basis of theory. This method applies 
not only to the normal conditions of flight, but applies 
equally well for flight with one or more engines inopera- 
tive. In case the power available is unsymmetrically 
placed because of failure of one or more engines, the 
power required is so greatly affected by the axes of 
thrust that the determination of performance of the air- 
plane for all conditions of weight, power, altitude, and 
temperature would be an endless process of empirically 
collecting the data, were it not for the correlation of the 
data with and by the means of theory. This method 
then is very excellently adapted to the flight test per- 
formance of airplanes within the present normal speeds. 
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The Column Strength of Closed, Thin- Walled 
Sections of 18-8 Stainless Steel’ 


HOWARD W. BARLOW 
New York Unwversityt 


SUMMARY 

Thin-walled members having closed cross-sectional shapes are 
used as columns in aircraft construction. The purpose of this 
investigation was to determine and evaluate the factors influenc- 
ing the behavior of such columns fabricated from 18-8 corrosion- 
resistant steel. Approximately 300 specimens comprising two 
different cross-sectional shapes in nominal thicknesses of 0.006 in., 
0.008 in., 0.010 in., 0.012 in., 0.014 in. were tested as columns 
under flat-end and under pin-end conditions. 

In order to carry out this investigation it was necessary to 
develop new methods which were found to have broad applica- 
tion. 

Two new types of special fixtures were developed which per- 
mitted pin-end tests to be conducted with the use of flat-end 
specimens. 

A simplified and improved technique was developed for testing 
the columns. This technique is superior to that currently em- 
ployed which makes use of imbedding plates, of cast ends and of 
specimens subject to loading in prior test. The technique in- 
cludes an improvement in the process of preparing the samples for 
test. 

The specimens investigated were obtained from the stock of a 
commercial manufacturer and had been made from two flanged 
channels by the shot-welding process. The channels, in turn, 
had been produced by the draw-bench method and were of full- 
hard, 18-8 chrome-nickel, corrosion-resistant steel. This ma- 
terial is austenitic at room temperature and cold working in- 
creases its strength, hardness and modulus of elasticity, and 
raises its proportional limit. The fact that the physical properties 
of the material are subject to change with cold working makes it 
necessary to determine the basic strength characteristics of the 
material after all processing has been completed. 

The compressive elastic properties of the completed sections 
were determined by the use of the ‘“‘Plug’’ method.'® This method 
permits the determination of the compressive stress-strain re- 
lationship for any hollow, thin-walled section up to values of 
average unit stress considerably greater than those at which 
secondary failure usually occurs. 

The column tests were analyzed on the basis of the compressive 
elastic properties of the material. For the sections investigated 
herein it was established that in the equation for primary failure 
by translation, the tangent modulus in compression represents 
the effective modulus. 

The coefficient of end-fixity for flat-end column tests was found 
to be 4.0. 

The system of non-dimensional coordinates advocated by W. R. 
Osgood was used in plotting the data from the compression stress- 
strain tests and the data from the column tests. The results for 
all gages of both cross-sectional shapes were found to lie on 
identical curves. 

It was established that the critical unit stress for secondary 
failure is directly proportional to R and is inversely proportional 
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to the square of the flat-pitch ratio. R is defined as a section 
constant for any given shape. The value of R does not vary 
with the thickness of the material and may be determined for 
each section by compression tests of a few short specimens 


INTRODUCTION 


rr ORDER that the stainless steels used in aircraft con- 
struction may be able to compete with the lighter 
alloys on a strength-weight basis it is necessary to re- 
duce the thickness of the individual sheets or sections 
to the minimum possible value. The fact that the 
elastic instability of these thin sheets varies inversely 
as the square of the flat-pitch ratio requires that flat 
surfaces be broken up by corrugations, ridges, vee 
grooves and other irregularities. Flat-pitch ratio is 
defined as the ratio of the width of the widest flat side 
to the thickness of the material. 

The practice of corrugating has led to the develop- 
ment of a large number of different cross-sectional 
shapes for stiffeners used in stressed-skin construction, 
and for columns used in truss-type construction. Some 
attempts have been made to investigate the column 
properties of such irregular cross-sectional shapes, but 
without a great deal of success. The hardness and 
toughness of the material, combined with the fact that 
it is used in such thin sheets, have presented difficult 
problems to investigators. 

When used as columns either as open sections, such 
as channels, zees, angles or hat-shaped profiles, or when 
used as closed sections of irregular cross-sectional 
shapes, such as are investigated herein, corrosion- 
resistant steel exhibits the same general traits as the 
lighter alloys, and is subject to the same general type 
of failure. However, the limitations accompanying 
these types of failure for the stainless steel have been 
subject to much less investigation than for the lighter 
alloys. 

The purpose of this investigation was to determine 
the fundamental principles underlying the column 
strength of hollow, thin-walled structural sections of 
stainless steel, such as are used, or may be used in air- 
plane or similar structures requiring extreme refine- 
ment of design in the matter of weight reduction, cor- 
rosion resistance and high strength. 


SCOPE 


The scope of this investigation includes column 
tests of two closed sections of typical shape made 
from sheets of nominal thickness 0.006 in., 0.008 in., 
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0.010 in., 0.012 in., 0.014 in. The material was 18-8 
chrome-nickel, corrosion-resistant steel. It is austenitic 
at room temperatures. Pin-end as well as flat-end tests 
were run. Compressive tests were made to determine 
the elastic properties of the material in compression. 
The scope was necessarily limited because of the limited 
number of shapes available for test. Thin stainless 
steel is not a material which may be readily worked by 
an amateur with the aid of a few hand tools. It is 
extremely hard and tough, and may be fabricated suc- 
cessfully only with the aid of special mechanical devices 
and by experienced workmen. 

Although the range of the investigation was limited 
to closed sections, the technique developed and used 
is also applicable to the testing of open sections such 
as mentioned above. Open sections are subject to the 
same type of failure as closed ones and in addition they 
may fail by twisting, as has been investigated elsewhere. ' 
A few tests were conducted using channels, in order to 
show that the special technique of preparation and 
testing could be extended to open sections. 

The cross-sectional shapes of the sections investigated 
are shown in Fig. 1. The section designations 5-22 and 
S-23 are those of Fleetwings, Inc., who furnished the 
material. The geometrical properties are given in 
Table 1. 





S-23 


Stainless steel sections tested. 


TABLE | 


Section Properties 


K x-x hey 





I y-y 


Sec. A I X-% 
S-22 2.228¢ 0.0614¢ 0.0722t 0.166 0.180 
$-23 3. 56 0.342t 0.310 


0.3245t 


0.302 


The areas and moments of inertia were calculated 
from the centerline dimensions of the sections. The 
moments of inertia were computed on the assumption 
that the corners were square, and the areas were com- 
puted on the same assumption, except that one thick- 
ness of material was deducted for each 90 degree bend. 
This latter procedure is consistent with shop practice, 
and the results obtained by it agree with those found 
by measuring flattened-cut specimens. The sections 
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are fabricated by drawing strips of full-hard, 18-8 
corrosion-resistant steel through a set of forming rolls 
which shape the flat stock into the desired contours. 
In order to produce closed sections, flanged channels 
are made. These are then joined at the flanges by a 
continuous roll-welding process which spaces the 
welds at any desired interval. The spacing for the 
sections used was approximately #/j, in. The speci- 
mens are uniform in shape and have the minimum 
possible amount of material along the joined flanges. 

The S-22 sections were furnished in nominal gages of 
0.006 in., 0.008 in., 0.010 in., 0.012in. The S-23 sections 
were furnished in nominal gages of 0.006 in., 0.008 in., 
0.010 in., 0.012 in., 0.014 in. The thickness of the 
actual material differed somewhat from the nominal 
value and was determined for each specimen by 
weighing. Owing to the compact shape of the cross- 
section it was difficult to determine the actual gage by 
measuring with a micrometer. The density of the 
material was determined from a piece of flat sheet, and 
found to be 0.290 Ib. per cu.in. 

It was not possible to run complete sets of tests on 
all of the gages listed above, since in some cases the 
supply of individual thicknesses was limited. This 
shortage did not hamper the investigation, however, 
and there are no appreciable gaps because of it. 

Prior to the undertaking of any actual column tests 
it was necessary to make a thorough study of the prob- 
lems involved in testing and in preparing for test the 
hard, thin-walled sections which were used. There is 
very little information on the subject, and that which 
was found was not applicable to the solution of the 
problems in hand. 


END FINISHING 


All of the specimens which were tested had flat, 
parallel ends. The matter of obtaining consistently 
similar ends within close limits required an improve- 
ment in the procedure currently used for this operation. 
It has been general practice to require that stainless 
steel column test specimens be finished to within 0.002 
in. of being flat and parallel, and the plane of each end 
be perpendicular to the longitudinal axis of the speci- 
men.* The conditions are usually obtained by grind- 
ing in a surface grinder. The use of a grinding wheel on 
the end of a clamped-down specimen is satisfactory in 
the case of heavy, thick-walled specimens, but is not 
applicable in the case of material 0.006 in. thick, for 
example. The thin elements vibrate under the action 
of the grinding wheel and a wavy surface is produced 
on the end of the piece. A great improvement in the 
method of finishing was obtained by grinding the ends 
of each specimen against an emery cloth disc which was 
fastened by rubber cement to a steel disc bolted to a 
lathe face plate. The specimen was held by hand in an 
angle bracket mounted in the lathe cross carriage. All 
parts of the apparatus were carefully machine-ground. 
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Fic. 2. End finishing apparatus for long specimens. 


The angle bracket was carefully set prior to the 
grinding of any group of specimens. The arrange- 
ment shown in Fig. 2 consists of a single disc, and was 
used for the longer specimens. In finishing any speci- 
men, one hand was used to hold it securely against the 
angle rest, and at the same time to push it against the 
emery disc. The other hand cranked the cross carriage 
of the lathe, carrying the specimen, back and forth 
across the semi-diameter of the disc in order to wear 
the disc uniformly. After two or three passes the 
piece was removed and the burrs taken off by a scraping 
tool. A finishing grind was accomplished by holding 
the specimen in the angle bracket as before, but allow- 
ing it to rub against the grinding disc only slightly. 
The cross carriage was again moved laterally in order 
to assure uniform wear. The piece was then rotated 
end-for-ehd about a horizontal axis and the process re- 
peated for the other end. In all cases roughing and 
finishing cuts were alternated. This resulted in keep- 
ing the surface of the wheel true, and permitted the use 
of one piece of emery cloth for over two hundred speci- 
mens. A careful operator can finish both ends of a 
specimen in about five minutes. 

The conditions of flat ends were checked by holding 
the specimens against a ground and hardened steel 
block held up to the light. If any light could be seen 
between the specimens and block, or if there was any 
tendency to rock, the specimen was refinished. A 
further check was provided for the short specimens as 
follows: the flat loading blocks used in the testing 
machines were carefully located, and shimmed up in 
order to make them parallel. The short specimens 
were, then inserted in the machine and the heads 
brought together in order to pinch the piece very 
lightly. If the specimen had flat, parallel ends normal 
to its longitudinal axis, it would rub smoothly at both 
ends when rotated by hand. Any irregularities were 
indicated instantly by the tendency of the piece to 
ride or rotate about any high points, or to resist rotating 
at all in the case of the ends not being parallel. 
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Owing to the fact that the ends of the sections tend 
to spring out of shape slightly, it was necessary to de- 
vise a two-disc grinding apparatus for the short speci- 
mens. When these pieces are flipped end-for-end, any 
minor distortion results in the production of ends which 
are not square and parallel after grinding. To correct 
this defect, a pair of the finishing discs was mounted on a 
ground shaft, and each disc was carefully and accurately 
ground. The specimen to be finished is held in a short 
angle bracket which, in turn, slides in the angle bracket 
of the original fixture. The use of the parallel grinding 
wheels removes the necessity of rotating the specimens 
in order to finish both ends. The arrangement is shown 
in Fig. 3. 
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Fic. 3. Two-disc end finishing apparatus for short 
specimens. 


CoLuMN Test DEVICES 


Pin-end column tests were conducted by means of a 
new type of column testing fixture which allows pin-end 
conditions to be obtained for specimens having flat 
ends.* This is shown in Fig. 4. The original device 
is made in pairs. Each unit consists essentially of a 
round, hardened steel bar supported with axis hori- 
zontal, between two ball bearing assemblies. The 
middle portion of the bar has a wide slot milled at its 
center. The surface of this slot is ground down to a 
plane which lies exactly in the plane of the diameters 
of the ball bearing assemblies. The non-Brinell rating 
of the fixture is 24,000 Ibs. and the flexural rigidity 
of the loading bar is such as to allow an angular deflec- 
tion of 1.2° under this load. The design play in the 
bearings allows an angular deflection of 1.5° in the 
loading bar prior to binding of the bearings. Since 
maximum test loads were less than 8000 Ibs. the 
fixture had an ample factor of safety against bearing 
failure or binding. 

The preliminary jig just described has many ad- 
vantages over other devices previously used for the 
same purpose. It is easy to install and align. The 
axes of rotation of the ends of the specimens lie in the 
plane of the loading and thus no correction for effective 
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Pin end column testing fixture (single degree of 
freedom). 


Fic. 4. 


length of the column is required. The loading bars 
may be statically balanced in the case of delicate 
tests, and pin-end tests may be run on flat-end speci- 
mens. 

Extensive use of the fixture’ indicated that it was 
quite satisfactory for the use for which it was designed. 





Pin end column testing fixture (two degrees of 
freedom). 


Fic. 5. 
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Furthermore, the scope of its application could be 
broadened if it were provided with two degrees of free- 
dom, instead of one, and redesigned so as to accommo- 
date specimens having a length less than six inches. 


A new jig, shown in Fig. 5, was built. This fixture 
has two degrees of freedom, uses half bearings, and can 
be used for the testing of specimens as short as one-half 
inch. It is counter-weighted in order to be statically 
balanced about both axes of rotation. Cradling pins 
are provided in order to prevent the individual units 
from falling apart during transit, during assembly in the 
testing machine, or when column failure occurs. The 
principles employed in this new fixture are basically the 
same as those used in the simpler one. 


Owing to the type of construction employed in the 
single degree freedom jig, it was not possible to conduct 
tests to determine the coefficient of friction, although 
the desirability of doing so was realized. Since the two 
fixtures were of similar construction it was believed that 
friction tests on the new type of fixture would be in- 
dicative of the results which might be expected from 
the simpler device. By removing the cradling arms 
and aligning plates it is possible to bring the two 
halves of the jig together for the purpose of measuring 
the coefficient of friction of the assembly. Owing to 
the extreme freedom of movement of the assembly this 
proved to be difficult in practice although it sounded 
well in theory. Since the bearing assemblies used for 
each axis of freedom were of the same size and type it 
was felt that a test about one axis would be representa- 
tive. This test was conducted by assembling the base 
castings and the accompanying half-bearings around a 
cylindrical steel bar. The torque required to start 
rotation of the bar under various compressiye loads 
on the assembly was determined. The coefficient of 
friction equals the torque divided by the product of 
the normal load multiplied by the diameter of ball 
centers. The effect of the friction is to produce an 
eccentricity which is equal to the coefficient of friction 
multiplied by the diameter of the ball centers. The 
results of the tests indicated that this would amount 
to an eccentricity of not more than 0.0039 in. at 100 
Ibs. and 0.00595 in. at 6000 Ibs. for the jig having a 
single degree of freedom. Later experience indicated 
that these were of negligible magnitude. In addition 
to these friction tests, two columns of '/2 in. diameter 
bars of cold-rolled steel 24 in. long, gave critical loads 
of 1523 and 1524 lIbs., respectively, when tested with 
the use of the single degree freedom device. On the 
basis of the accepted value of the modulus of elasticity 
of 29 X 10° Ibs. per sq.in., this would indicate an end- 
fixity coefficient of 1.0 for the fixture. 


From the results of the two types of tests discussed 
above, it must be concluded that the fixtures described 
and used in the investigation of the stainless steel 
columns would provide as close an approach to the 
ideal value of C = 1.0 for pin-end tests as is necessary 
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for accurate experimental testing and development 
work. 

The pin-end column testing fixture described above 
represents one of the best devices now in use for testing 
columns. Notwithstanding the existence of other 
devices developed for the same purpose, various investi- 
gators have preferred to conduct column tests with 
flat blocks. It has been assumed that flat-end tests 
develop end-fixity coefficients having values ranging 
from two to four, and in some instances it has been 
assumed that the restraint is the same as that developed 
for contiuuous members in aircraft structures. 

The lack of definite information on the end-fixity 
produced by flat-end tests seemed to be a good reason 
to investigate and evaluate it with the aid of the pin- 
end testing equipment. 

A complete series of flat-end column tests was run on 
both the S-22 and the S-23 sections; these were dupli- 
cates of the pin-end tests. Fig. 6 shows a specimen 
under test. 

The column tests were conducted on two different 
Tinius Olsen, hand-operated, screw-type machines. 
One was located at New York University and used for 
testing columns ten inches or less in length, while the 
other machine was located in the structures laboratory 
of Fleetwings, Inc., and used for testing the longer 
columns. 


IMBEDDING PLATES 


It has been common practice to use imbedding plates 
between the ends of stainless steel column test speci- 
mens and the steel loading blocks of the testing ma- 
chine. Lead, copper, aluminum, aluminum alloy, 
wet cardboard, pyralin sheets softened in hot water, 
and various other compositions have been used. The 
purpose of the plates is to restrain the ends of the test 
specimen and thus to avoid failure or distortion at the 
ends. It has been assumed that this practice would 
produce a uniform stress distribution over the ends of 
the specimens. 

Several tests with aluminum alloy imbedding plates, 
as well as without imbedding plates were conducted on 
duplicate specimens. Flat-end test procedure was 
used. It was found that there was no local failure at 
the ends when no plates were used, but that the use of 
imbedding plates gave unsatisfactory results. In the 
latter cases the thin outstanding flanges of the specimens 
resisted imbedding and buckled locally. Furthermore, 
the results from the tests with imbedding plates were 
erratic, as compared with those from tests without the 
imbedding plates. This experience indicated that the 
use of imbedding plates introduced another cause for 
failure and was quite unnecessary. 

The only satisfactory method for securing even 
distribution of load is to have the ends finished flat and 
parallel. Neither imbedding plates nor any other de- 
vices will remedy the faults of improperly finished ends. 


THIN-WALLED SECTIONS 155 





Typical flat end column text (long column). 


CASTING THE ENDS 


It has been suggested by some investigators that the 
ends of the thin metal sections might be cast in small 
slugs of material such as Wood’s metal (melting point 
about 150°F.). The purpose of this is to provide a 
restraint for the ends of the specimens and prevent 
any local distortion. In order to check the value of this 
procedure, a number of samples of open and of closed 
sections were prepared with and without cast ends, 
and flat-end tests were conducted with them. The 
results of the tests on the closed sections were generally 
erratic. The results of flat-end tests on a series of 
1/,in. X '/2in. X 0.010 in. channels indicated that the 
major effect of casting the ends of this section is a 
shortening of the effective length of the column by an 
amount equal to the portion of the column encased in 
the casting. Casting made no difference in the manner 
of failure. It may be concluded, therefore, that this 
process is not necessary, and in fact, is likely to indicate 
a false value of the column strength. 


USE OF PRE-STRESSED SPECIMENS 


In some cases it has been the practice to test a piece 
having a high slenderness ratio and, after failure, to 
cut two shorter specimens from the apparently un- 
damaged portion of the column. This is done in order 
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to obtain three different tests from the one long speci- 
men. The impropriety of this as applied to structural 
steel columns has been noted by Basquin.® In order to 
check the inapplicability of this procedure to the stain- 
less steel testing, several flat-end tests were conducted 
using these so-called pre-stressed specimens. These 
had all been cut from longer pieces previously tested. 
In most cases it was found that the pre-stressed speci- 
mens of a given length failed at higher loads than new, 
unstressed specimens of the same length. As in the 
case of the specimens tested with imbedding plates, 
the results from the pre-stressed pieces were generally 
erratic. The conclusion was confirmed that it was not 
advisable to use specimens which had been subject to 
loads in previous tests. 


CENTERING UNDER LOAD 


Current practice in column testing makes use of 
centering the compressive specimens under load.®*? 
The purpose is to account for any irregularities or devia- 
tions from straightness of the column, by means of 
centering the column and adjusting the location of its 
ends. These are moved until the lateral deflection of 
the middle of the column is extremely small and does 
not vary between loads of 150 Ibs. and loads represent- 
ing about 90 per cent of the maximum load expected. 
The refinement of the column testing process as repre- 
sented by this method increases the accuracy of the 
results, but in some cases may be superfluous. 

The properties of the stainless steel used in this in- 
vestigation were such as to preclude the use of this 
procedure. It is inadvisable to impose any degree of 
pre-stressing on this material because of its extreme 
sensitivity to cold working. This response to cold 
working is generally advantageous, since it increases the 
strength and elastic properties of the material. How- 
ever, varying amounts of cold working have varying 
effects, and the repeated loading to high stresses in- 
volved in the process of centering under load would 
undoubtedly result in a greater scattering of test points 
than would the small eccentricities obtained in original 
centering, and the small deviations from straightness. 
Moreover, Osgood, a proponent of centering under 
load,’ says, “If, when the free lengths of the test 
columns have been computed, plots of average stress at 
failure against ratio of slenderness based on the com- 
puted free length, or modifications of these quantities 
that take into account the variations in the properties 
of the material, lie on a reasonably smooth curve, this 
fact may be regarded as proof that small initial curva- 
tures under conditions of proper centering do not affect 
the strengths of even inelastic columns freely supported 
or elastically restrained at the ends.’’ The fact that 
such consistent results were obtained in this investiga- 
tion substantiates the soundness of the centering process 
which was used. 

For the flat-end tests this procedure consisted 
simply of locating the loading blocks approximately in 
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the center of the testing machine heads, and locating 
the specimen, by eye, in the approximate center of the 
loading blocks. For the pin-end tests, the loading bars 
were carefully lined up in the machine so that their 
axes were in the plane normal to the lower head of the 
testing machine. The specimens were centered with 
respect to the bars by a careful procedure which re- 
stricted the error to 0.0005 in. 


LATERAL MOVEMENT OF MOVABLE HEAD 


Basquin,® in his study of tests on structural steel 
columns, noted that there might be considerable dis- 
crepancy due to tilting or shifting of the heads of the 
testing machine. Means were not available for pre- 
venting these motions on the machines used here, but an 
effort was made to have the specimens centered in the 
machine. It is not believed that these variables were 
of much importance in the tests reported in this paper. 
The test points show such a small tendency to scatter 
that the effect of unrepeatable factors, such as testing 
machine head tilting or shifting, must have been 
negligible. 


TEST PROCEDURE 


As far as possible in view of the supply of the various 
gages of the two sections, S-22 and S-23, pin-end tests 
were conducted on specimens having slenderness ratios 
from 20 to 140, in intervals of ten. A slenderness ratio 
of 140 for the S-23 section represents a length of 42 
in. A duplicate set of specimens having a similar range 
of slenderness ratios was tested as flat-end columns; 
in addition, a number of shorter pieces, down to lengths 
of one inch, were tested as flat-end columns. All loads 
were applied slowly and evenly. The tests were con- 
ducted on hand-operated machines; the beam was 
always kept balanced, and the critical load determined 
by the drop of the beam. 


ELASTIC PROPERTIES OF THE MATERIALS 


‘“‘Plug’’ method tests were conducted in order to de- 
termine the elastic properties of the material in com- 
pression.!5 Nominal gages of 0.006 in., 0.008 in., 
0.010 in., 0.014 in. were tested for the S-23 sections and 
of 0.006 in., 0.008 in., 0.010 in., 0.012 in. were tested 
for the S-22 sections. Huggenberger type “‘A’’ tensome- 
ters were used to measure strain. Two instruments 
were used for each test. They were attached near the 
midpoint of the sections at points near the ‘‘outside”’ 
corners of the flat sides having the greatest width. 
Loads were applied slowly and evenly, and frequent 
readings of strain were made. The movement of the 
machine head was stopped each time a strain reading 
was taken, and the indicating hand of the strain gage 
was allowed to come to rest before reading. The process 
permitted a condition of equilibrium to be established 
in the piece. 

Duplicate specimens of each gage and section were 
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tested. In about one-half of the tests the loading 
process was stopped at a point representing approxi- 
mately the yield strength of the material, as defined 
later. At this point the load was allowed to decrease. 
Strain readings at various loads were taken during this 
unloading process which was continued until all load 
had been removed. In the remainder of the tests the 
load was applied until secondary failure occurred, 
as indicated by a wrinkling of the free flange or of the 
widest flat side. The results of the material tests are 
shown in Fig. 9. 


COLUMN THEORY 


Lundquist proposes a rational classification of column 
failure which is dependent upon the type of failure which 
occurs."* In this classification, column failures are 
divided into two groups: 

(a) Primary, or general failure, which is defined as 
any type of failure in which the cross-sections are 
translated, rotated or both translated and rotated, 
but not distorted in their own plane. 

(b) Secondary, or local failure, which is defined as 
any type of failure in which cross-sections are dis- 
torted in their own plane, but not translated or rotated. 

It was found that analysis of the columns investigated 
herein could be based upon this classification of Lund- 
quist’s. The phenomena related to these types of 
failures are generally well known. For primary failure, 
the equation for critical stress is: 


ad 
P/A = cr 


_ =. 1) 
(L/Ky* ' 





where £’ is used to designate the effective modulus of 
elasticity. It is customary to write this equation as 
follows: 





tA = a *E" = 2) 
P/A (L’/K)? (2, 
where L’ equal L/+/c. 

For cases where the critical stress for primary failure 
is less than the proportional limit of the material, Eq. 
(2) is the Euler equation,’ and E’ has the value of £,, 
the Young’s modulus in compression. For columns 
where the stress at primary failure exceeds the propor- 
tional limit of the material, Euler’s formula no longer 
applies. For such cases the ‘double modulus’ is 
applicable.® The double modulus has a value some- 
where between the Young’s modulus in compression 
and the tangent modulus in compression, £,, at the 
particular critical stress. The exact value depends 
upon the material as well as upon the cross-sectional 
shape of the piece. Templin established that the 
tangent modulus could be used in place of the double 
modulus for various aluminum alloys,'' and Basquin 
reached the same conclusions after analysis of steel 
column tests.6 This procedure also proved to be valid 
for the sections and materials investigated herein. 
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Since the end-fixity coefficient, C, for the pin-end 
tests is known to have a value of 1.0, it was possible 
to determine the value of the end-fixity coefficient for 
the flat-end test specimens from the flat-end test results. 
It was found that in all cases, for both sections and all 
gages the value of this was four. 


Non- DIMENSIONAL COORDINATES 


Whenever column test results are represented by 
curves wherein average stress is plotted against slender- 
ness ratio there will be some unavoidable scattering of 
the test points. This variation is due to the variations 
in the geometry of the individual sections as well as 
the variations in the material from which they are made. 
This latter variability can be corrected to the extent 
to which the stress-strain relations remain invariable 
along the length in any one given specimen and are 
affinely related from specimen to specimen.’ Osgood 
proposes that the stress-strain relationship be reduced 
to a non-dimensional relationship,’ and also proposes 
that column test results be likewise plotted on non- 
dimensional coordinates. 

Considering this latter suggestion first, it is noted 
that Osgood proposes the following coordinates :’ 


og, = (P/A)/S as ordinate, and (3) 


d (1/m)(L'/K) VS/E, as abscissa (4) 


‘os 


It was suggested in reference 7 that S be determined 
from the compressive stress-strain diagram of the 
material in each specimen as the intersection with the 
stress-strain curve of a line through the origin having a 
slope BE where £ is a constant for a given material, and 
varies between 0 and 1.0, and E is taken as £#,, the 
Young’s modulus in compression. 

The value of S may be chosen also in ways different 
from that proposed above. In this paper S is taken as 
the compressive stress at which the material exhibits a 
shortening under load of 0.002 in. per inch in excess of 
that which would be computed from a Young’s modulus 
of elasticity, Z, equal to 25 X 10°, and the usual 
formula: o = E, X «. 

The value of the stress S so obtained corresponds to 
the definition of yield strength given in Navy Depart- 
ment specification 44, T 27, and is sometimes referred 
to as proof stress.’ If this value of S is used, the 
corresponding value of E, of 25 X 10® should also be 
used in the non-dimensional relationship. 

In some cases it may be more convenient to select S 
arbitrarily as a value close to the yield strength of the 
material (as previously defined). This latter method 
will generally give satisfactory results provided the 
value of S so selected does not vary greatly from speci- 
men to specimen and represents a rough mean of the 
yield strengths. Actually it does not make much dif- 
ference how S and £ are selected, provided E has a 
value close to the compressive modulus of elasticity of 
the material at zero stress, E,, and provided S has a 
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value approximately equal to the compressive yield 
strength of the material. 

The advantage of plotting on the non-dimensional 
coordinates is that all column test data or all stress- 
strain data, whichever is being plotted, will fall on one 
and the same curve, provided the various original] 
stress-strain curves of the material are affinely related. 

It follows as a logical conclusion that these individual 
curves, one developed from the material test results 
and the other developed from the column test results, 
will be identical provided the proper values of the 
variables and the constants are selected. It is neces- 
sary, of course, that the same system of non-dimen- 
sional coordinates be used for plotting column test 
results as well as material test results. 

In reducing the column test results for primary failure 
to the non-dimensional system, reference is made to 
Eq. (3), where P/A should be taken as the critical 
average unit stress at which primary failure occurs. 

In order to reduce the material test data to the non- 
dimensional system, Eq. (4) must be modified in such 
a manner that it will include the compressive elastic 
properties of the material. If the value of (L’/K), as 
obtained from Eq. (2), is substituted in Eq. (4) the 
following relationship is obtained: 


(E) _ (S) Ys 
Mas = (1/xr)(x x " = ia x (5) 
) (P/A) (E.) P/A E 
In order to correlate column test data with material 
test data, it is necessary to select the proper value of the 


effective modulus £’, which when substituted in Eq. (5) 
and used with Eq. (3) will produce a curve identical to 





Type of secondary failure—Sec. S-23. 
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that produced from the use of Eqs. (3) and (4). This 
value of E’ is obviously the value which satisfies Eq. 
(2), which is the general formula representing critical 
stress for primary column failure. In accordance with 
the currently accepted theory, E’ represents the double 
modulus, E, which for some materials and sections, 
including those investigated herein, may be in turn 
represented by the tangent modulus in compression, 


E,. 


SECONDARY FAILURE 


For the analysis of secondary failure Timoshenko’s 
methods may be applied.’ The phenomenon of second- 
ary failure of sections having flat elements depends 
upon the elastic stability of the individual elements 
composing the section, as well as upon the degree of 
edge support furnished to each element by all or any 
of the other elements of the section. 

The critical stress for secondary failure is defined 
as the average unit stress over the whole cross-section 
existing when any element of the section becomes elas- 
tically unstable. If the load is increased after secondary 
failure occurs, the stress distribution over the section 
will no longer be uniform and with further increase in 
load additional elements will fail until the entire section 
collapses. This progressive failure is probably some 
function of the effective slenderness ratios of the various 
elements of the cross-section. 

The critical stress for elastic stability of a flat plate 
subject to end compression and having various degrees 
of support along the edges is given as: 


kr°t?E, 


wane (6) 
12(1 — p?)b? 


Sav. - 
This equation as written applies only to stresses below 
the proportional limit of the material. It appears logical 
to apply it to conditions existing beyond the propor- 
tional limit by substituting the effective modulus in 
compression E’, for Young’s modulus E,. The equation 


is then: 
kn*t?E’ : 
a 7) 
For 12(1 — y?)b? 


For sheets having a large ratio of length to width, k 
is a factor dependent only on the degree of restraint or 
elastic support provided at the edges of the sheet. 

For the purpose of analyzing entire sections of com- 
plex shape by investigation of the stability of the flat 
sides, the term k as used by Timoshenko for the simpler 
cases of flat sheets having various degrees of support 
along the edges may be replaced by a particular value 
R. This R is defined as an effectiveness factor. It 
represents the effectiveness of the restraint and sup- 
port provided at the edges of the element which fails 
first. It applies only in the cases where the length of 
the flat element is several times its width. 

It was determined in this investigation that R is sub- 
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stantially constant for a given geometrical section and 
invariant with thickness of material. 

The value of R may be determined from the crinkling 
tests of a few short specimens. From results of such 
tests and by means of Eq. (7), R is found from the 
following expression : 


_ (P/A)12(1 — y?)b? 


— (8) 
rk’? 


R 
The maximum slenderness ratio at which secondary 
failure is critical may be determined by equating the 
critical stresses for primary and secondary failure Eqs. 
(2) and (7) wherein R is substituted for k in Eq. (7), 
then, 
f mE’ = RE’? 
7 “" (L'/K)2— 12(1 — p?)b? 


from which: 
L'/Ke.=— (9) 


Typical forms of secondary failure are shown in Figs. 
7 and 8. 


TEST RESULTS 


Data from the compressive stress-strain tests of the 
fabricated sections as shown in Fig. 9 were reduced, 
by using Eqs. (3) and (5), to the non-dimensional basis 
previously referred to. Results of the calculations are 
plotted in Fig. 11 and a faired curve is drawn through 
the points. 





Fic. 8. Type of secondary failure—Sec. S-23. 
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Typical results of column tests are shown on Fig. 10 
where average unit stress is plotted against slenderness 
ratio. Data from tests of the other section and the 
other gages gave similar results. Flat-end tests have 
been reduced to equivalent pin-end lengths on the 
assumption that C for flat-end tests is 4.0. The curve 
drawn on Fig. 10 represents Eq. (2) wherein the tangent 
modulus in compression £,, as obtained by the ‘‘Plug”’ 
method, is used for the effective modulus E’. The 
agreement of the results of column tests with the curve 
so drawn is a first substantiation that the tangent 
modulus in compression so obtained satisfies the equa- 
tion for primary failure of the columns tested in this 
investigation. No attempt has been made to analyze 
the type of failure that occurred after secondary failure 
took place. The horizontal line on Fig. 10 represents 
the value of the average unit stress at secondary 
failure as obtained by the use of Eq. (9). This will be 
discussed later. 

The results of column tests on all gages of both sec- 
tions have been reduced to the non-dimensional values 
represented by Eqs. (3) and (4), where the values of S 
and £, for each gage and section are the same as those 
obtained in the tests of reference 15. Flat-end test re- 
sults are reduced to pin-end test results on the basis of 
C = 4.0 for flat-ends. The test points representing 
column test results are plotted on Fig. 12. The curve 
drawn on this figure is identical with that of Fig. 11. 

Values of o,, representing secondary failure in ac- 
cordance with Eq. (9), are shown as horizontal lines 
on Fig. 12. The value of } as used in this equation is 
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Fic. 10. Typical result of column tests 
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the width of the flat side which fails first. This width 
is measured from midpoint to midpoint of the adjacent 
supporting flanges (or the point of tangency of large 
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radii, as for the S-22 section). The value of » was as- 
sumed to be constant and equal to 0.30. The factor R 
was evaluated from tests of short columns and deter- 
mined to be as follows: 


R = 3.60 
= 6.90 


» 
| 


Critical values of slenderness ratio, L’/K, expressed 
by Eq. (9), were found for the nominal gages of the two 
sections. These, in turn, were substituted in the di- 
mensional and non-dimensional equations for column 
strength and plotted on the diagrams which show the 
column test results. In order to prevent confusion in 
the regions of the curves representing low values of 
As, many of the points representing tests of very short 
columns were omitted. There were so many of these 
that any plotting of them would only obscure the im- 
portant test points. 

There is presented in Table 2 a summary and com- 
parison of critical stresses for secondary failure as deter- 
mined from four different sources, as follows: 

1. From the column curves, such as shown on Fig. 
10, where P/A is plotted versus L’/K. The critical 
stress for secondary failure may be obtained by using 
the critical value of L’/K from Eq. (9) in conjunction 
with the column curve. 

2. From the non-dimensional 
drawn using particular values of S obtained from the 
compressive stress-strain curves for the particular sec- 
tion and gage. The critical stress for secondary 
failure is found from Fig. 12. o, and X,, are obtained 
by the use of Eqs. (3) and (4). The value of L’/K 
used in Eq. (4) is the value obtained from Eq. (9). 

3. From the non-dimensional curves which were 
drawn using an arbitrary value of S equal to 125,000 
Ibs. per sq.in.'® The critical stress for secondary failure 
is then found in the manner given for item (2) above. 

4. From the results of a few tests of short columns. 


curves which are 


CONCLUSIONS 


1. An inexpensive and rapid method has been de- 
veloped for preparing flat-end specimens of thin-walled 
stainless steel for column test. 

2. A fixture having a single degree of freedom of 
rotation has been developed for obtaining pin-end 
conditions in column testing. This device permits 
accurate reduction of column data and gives consistent 
results as judged from the close correlation obtained 
with material test data. A modification of this fixture 
has been developed which will permit obtaining two de- 
grees of freedom. Both fixtures permit the use of flat- 
end specimens and require no correction for effective 
length. 

3. Evidence indicates that the use of imbedding 
plates in column testing is unnecessary, provided the 
ends of the specimens have been accurately finished 
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TABLE 2 


Summary and Comparison of Critical Stress Values for Secondary 





Failure 
R—Cale. Act. 
1 o2 a3 O45 from Ga. from 
Sec. ————~+ 1000—————__ Test Test 
S-22-006 59.0 66.8 
008 92.5 88.2 93.7 
010 i18.0 114.0 114.0 
012 127.0 129.1 127.5 
( >44.0 >7.15 0.00641 
S-23-006 37.6 <52.8 <9.2 0.00641 
=37.8 ~6.75 0.00597 
>56.2 >6.30 0.00801 
= 24 23 OA <63.5 <8.04 0.00801 


010 81.5 81.7 80 
014 121.0 120.6 121.5 117.4 6.80 0.01369 





flat and parallel. The use of imbedding plates is likely 
to produce premature secondary failure. 

4. The use of pre-stressed specimens gives unsatis- 
factory results, and there is little justification for using 
material which has been subject to prior test. 

5. There is no need to cast the ends of column test 
specimens in Wood’s metal or other material in order to 
prevent secondary failure at the ends of the specimens. 
Adequate care in finishing the ends of the specimens 
has no substitute. 

6. The ‘Plug’ method has been developed for use 
in determining the compressive stress-strain relation- 
ship for thin-walled sections. It is especially applicable 
to built-up or other fabricated sections of material for 
which the elastic properties may have been altered by 
cold working or other processing. 

7. Carefully conducted column tests with flat-end 
specimens on flat, parallel loading blocks will develop 
end-fixity coefficients C, of 4.0. 

8. Thin-walled closed sections of 18-8 stainless steel, 
of the type investigated herein, develop primary 
failure (becomes elastically unstable as a whole) at 
stresses indicated by Eq. (2), where E’, the effective 
modulus of elasticity, is represented by the tangent 
modulus of elasticity Z;. This tangent modulus must 
be determined from compressive stress-strain tests on 
the fabricated sections, particularly if the material is 
such that the elastic properties are affected by cold 
working or other processing incident to forming. 

9. Secondary failure of the above type of specimens 
occurs at a critical effective slenderness ratio expressed 
by Eq. (9), in which R is a section constant for any 
given geometrical shape, and is invariant with the 
thickness of the sheet. 
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10. Results of column tests on sections such as are 
investigated herein may be successfully represented on 
plots using non-dimensional coordinates whose values 
are expressed in Eqs. (3) and (4). The results of these 
tests when so plotted show good correlation with the 
results of compressive stress-strain tests when ‘the 
latter are similarly reduced to the same basis with the 
aid of Eqs. (3) and (5). The compressive tangent 
modulus should be used to represent the effective modu- 
lus in all of the expressions for primary failure. 
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Shear Center of a Multi-Cell Metal Wing 


H. W. SIBERT 


Unwersity of Cincinnati 


SUMMARY 


Simple formulas based on the least work of shear are given for 
the coordinates of the shear center and the shear force per unit 
length of arc of a multi-cell wing beam. The terms involved in 
these formulas are ordinary summations or can be evaluated by 
Simpson’s Rule according as concentrated areas or skin effective 
in bending is concerned. A numerical example is worked out for 
the case of a two-cell wing beam with skin on the tension side 
effective in bending. 


INTRODUCTION AND NOMENCLATURE 


ban GENERAL principles for obtaining the shear 
center of a multi-cell metal wing by means of the 
least work of shear have been given by Hatcher.! Due 
to the complexity of the analysis only single-cell beams, 
such as leading-edge wing beams, have been treated so 
far.» ** 4 The purpose of this paper is to extend the 
methods developed in references 2 and 4 to metal wings 
containing two or more cells. Since the effect of taper 
on shear center is apparently small for modern metal 
wings,* the analysis will be limited to a wing of uniform 
cross-section. 

The following symbols and sign conventions will be 


used in this paper: 


cc = cross-section for which shear center is de- 
sired 

s = variable point on the skin of C, 

s = distance along the skin of C, from a starting 
point to point S (always +) 

Ys = shear flow, or shear force/unit length of skin, 
at S (+ when in direction of arrows shown 
in Fig. 2) 

I = moment of inertia of entire cross-section C, 


with respect to the principal axis about 
which the beam is bending 
L,t = length and thickness of a segment of the skin 
of C, 
Kun = J,y ds/t = SL/t for segment MN (always 
+) 
Qun = Statical moment about a principal axis of 
that part of the skin of C, between points 
M and N which is effective in bending 
= Si Qysds/t = Sy QustKus, M and N 
being points 
P = a fixed point on web ca of C, between points 
cand a (Fig. 2) 
Ayn = area bounded by the portion of the skin of 
C, between points M and Nand by the 
radial lines PM and PN (always +) 


J MN 
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Byun = Bing Qus dA ys, M and N being points 

a; = area of a concentrated area whose center of 
gravity is at a point z 

Mp = moment about P of the shear flow in the 
entire cross-section C, (+ when counter- 
clockwise) 

Dp = distance from P to the shear center meas- 


ured parallel to the principal axis about 
which the beam is bending; D = Mp/F 
Y, + Y, + 4>Y; (odd 2’s) + 2>-Y; 
(even 7’s), which is the summation in Simp- 
son’s Rule when Y is the integrand of the 
integral to be evaluated and the interval 
between the limits of integration is divided 
into (even) equal parts 
Subscripts 1, 2, 3 when applied to K, Q, and A refer to 
the closed circuits abca, acdea and dfged, respec- 
tively, of Fig. 2. 


Disr Y; = 


FUNDAMENTAL FORMULAS 


For constant cross-section and constant shear modu- 
lus the work of shear in the wing between C, and the 
loaded cross-section is‘ 


W = (B/2G) [qs* ds/t (1) 


where G is the shear modulus and B is the spanwise 
distance from C, to the loaded cross-section. 

Fig. 1 shows a portion of the beam of spanwise length 
dx with C, the face at distance x from some origin. 
With the shear flow in C, as shown in Fig. 1, the summa- 
tion of forces in the x direction gives gq, dx + gq dx = 
gsdx + qydx + dP, where P is the resultant of the bend- 
ing forces on the portion of C, shown in Fig. 1. Let F 
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Fic. 1. Relation between shear flow and bending forces. 
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be the component of total shear force on C, perpen- 
dicular to the principal axis about which the beam is 
bending and let F be positive when in the positive direc- 
tion of the other principal axis. Then,dP/dx = FOQ/I, 
where (Q is the statical moment of the portion of C, shown 
in Fig. 1. Hence, gi: + @ = qs +q1 + FQ/I. Note that 
gq: and gq can be thought of as flowing into, and g; and q 
as flowing out of, the portion of C, shown in Fig. 1. 
Then, for the particular case of F = 1, the last relation 
connecting 41, ge, gs, and g, can be written 


d¢ (in) = dog (out) + (1/J) 20 (interior) (2) 


Since g (in) and q (out) act along the cuts which bound 
the sub-section of C,, Eq. (2) is also valid when the sub- 
section of C, contains a closed contour of the skin. 
Thus, Eq. (2) is valid for all sub-sections of C,, provided 
that F = 1. 

In reference 4 it was shown that the contribution of a 
segment MN of C, to the moment Mp about a fixed 
point P equals the integral of (2g; dA ys) from M to N. 
Thus, Mp is the summation of these integrals taken 
over all the segments of C,. When Eq. (2) is used to 
find gs, Dp = Mp because Dp = Mp/F and Eq. (2) 
was obtained by taking F = 1. Hence, when gs is 
found from Eq. (2), 


Dp/2 = Mp/2 = S(S," Qs dAus) (3) 


in which the summation is taken over all segments 
MN with each term positive or negative according as 
the moment of the shear flow in that segment is counter- 
clockwise or clockwise, because Mp itself is positive 
when counterclockwise. 

Let x and y be the principal axes of C, and let x point 
to the right and y point upward. Since Dp = Mp/F, 
Dp is positive when both M>p and F are positive, that is, 
when Mp is counterclockwise and F is up for bending 
about the x axis or to the right for bending about the y 
axis. Hence, a positive Dp in Eq. (3) requires the 
shear center to be to the right of the y axis or below the 
x axis according as there is bending about the x or y 
axes. 


THREE-CELL WING 


The three-cell beam in Fig. 2 is triply indeterminate 
so far as shear flows are concerned, and the shear flow gs 
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Fic. 2. Three-cell wing. 


at any point can be obtained only after the shear flow 
at one point in each of the three cells has been found. 
Let the three shear flows to be determined be q,, g,’ and 
gq at infinitesimal distances to the left of point c, to the 
right of c, and to the right of point d, respectively (Fig. 
2). By means of Eq. (2) and suitable sub-sections of 
the beam it can be shown that gs has the following 
values: 


ds = & — Qs/I on cha 

Is = & — & + Qs/I on ca 
ds = a’ + Qs/I on cd 

Is = a’ — Ga + Qeas/I on de 
ds = ga + Qas/I on dfge 

Qs = &’ — (Qi + Qas)/I on ae _] 


(4) 





The three least-work conditions, O0OW/dg, = 0, 
OW/0q.’ = 0 and OW/0qz = 0 can now be applied to 
Eq. (1). By virtue of Eq. (4) and the definitions of 
K and J these three conditions yield, respectively, 


Kige — Kacge’ = (Seta — Sca)/I (5) 
Keg.’ — Kacge — Kaede = 

(QKae + Seca + Jace — Scée)/I (6) 
Ksqa — Kaege’ = (QeaK ae + Jae — Sayge)/I (7) 


the term Q,4K4, in Eq. (7) being due to the breaking 
up of Q.g5 on segment de into Q.4 + Qgs. Note that 
Yer Qc’, and gy can be obtained by solving Eqs. (5), (6), 
and (7) simultaneously. 

From Eggs. (3) and (4) and the definitions of A and 8, 


Dp/2 = Aige + A2g.’ + Asga + 
(Bede + Barge ae Biba vies Bae —_ Q:A,.)/1 (8) 


because with point P on ca the moment of the shear 
flow in ca about point P is zero and A,,, = Ay, Acg + 
Ag + Aa = Ag and A aye = Aw = A. Eq. (8) 
gives the correct sign of Dp because the values of gs in 
Eq. (4) were based on the directions of shear flow 
shown in Fig. 2 and these directions all give a counter- 
clockwise moment about ‘point P when P is on ca. 

All the formulas necessary for locating the shear 
center of a three-cell beam have now been given. A 
similar analysis will solve an n-celled wing (m > 3). 
For such a wing there will be m simultaneous linear 
equations connecting the unknown shear flows, g,, 
Gc’, Yq -.-. in the n cells. 


Two-CELL WING 


Let the two-cell wing be the wing in Fig. 2 with sec- 
tion dfge cut off. The effect of removing section dfge is 
to eliminate entirely Eq. (7), which came from 0W/0q, 
= 0, and to delete from Eqs. (4), (5), (6), and (8) all 
terms containing g, or pertaining to section dfge. The 
final results are: 
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— Jea)/I (9) 
in J cae) /1 (10) 


Kiq = Kaede’ - (Seba 
Keg.’ im Kacd are (QKae + Jes + Jes 


Dp/2 = Aiq + Ag,’ +}- 


(Bede ‘ie Biba ey Bae a Q:Aq-)/I (11) 


After g, and g,’ have been found from Eqs. (9) and (10), 
Dp can be obtained from Eq. (11). 

An alternative form for Eqs. (10) and (11) can be 
derived from the first two formulas of Eq. (4) and the 
relation gs = 4.’ + Q.s/IJ on cdea. Then, the two 
least-work conditions, oW/0q, = 0 and oW/0q., = 0, 
applied to Eqs. (1) and (3) change the right-hand side of 
Eq. (10) to (Iza — Jedea)/I and the last term of Eq. 
(11) to (Bigeo — Bera)/I. By means of the relations 
Qi + Qidea = 0, Jun + Jnu = QunKun and Byyw 
+ Byu = QunwAwy in conjunction with Eq. (17) and 
the corresponding equation for 8gcpz these alternative 
forms of Eqs. (10) and (11) can be transformed into 
the original forms. 

By definition the first power only of Q is involved in J 
and 8. Moreover, the Q for the skin alone is zero for 
any portion of the skin not effective in bending. Hence, 
the total Q, J or 8 for any segment of C, equals the sum 
of the Q, J or 8 of the concentrated areas associated 
with the segment plus the Q, J or 8 of that portion of the 
skin of the segment which is effective in bending. 


SKIN EFFECTIVE IN BENDING 


In reference 4 formulas have been given for J and 6 
for a segment of the skin of constant thickness when it is 
effective in bending. If MN is such a segment and if 
the segment is divided into m (even) equal parts with the 
points of division numbered from 0 to n beginning with 
point M and ending with point N, the 8, J and Q for 
bending about the x axis are: 


Bun = Lyn tun(dosrAwiyi)/3n (12) 


Jun = {Cy-2 + 2C,-4 + 3C,-6+... 
+ (n/2)C)}4L?%yy/3n? (13) 
Qun = Lun tun(Disryi)/3n = 
(Cot Cr+ Cyt... +C,)2Luntun/3n (14) 
where 
Co = 1 + ¥0/2, Cy = Yu—1 + In /2 
(15) 


C; = ¥i-1 +H + N41 A Norn 


If a J (or 8) refers to a segment comprising two or 
more sub-segments, this J (or 8) can be reduced to the 
sum of the J’s (or 8’s) on the sub-segments plus certain 
quantities involving Q. This is accomplished by letting 
the Q on each sub-segment equal the sum of the Q’s on 
the previous sub-segments plus the Q on the segment in 
question; for example, Qgcs on sub-segment CD of a 
segment BCD equals Qgc + Qcs. By virtue of such 
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relations and the definitions of K and J it can be proved 
that 
(16) 


(17) 


J Bcp 7a J Bc . Joo a QOscKep 


Jacpor = Jscv + Joe + OscoKpe 


The corresponding formulas for Bgcp and Szgcpex are 
obtained by replacing each J and K in Eqs. (16) and 
(17) by 6 and A, respectively. By means of these rela- 
tions the problem of finding any J or B in Eggs. (5) 
through (11) reduces to that of finding the J or 8 on 
segments of constant thickness, to which Eqs. (12) and 
(13) are applicable. Note that when Q is zero at every 
point of a segment, J is also zero for that segment. 
Thus, Jgcp = Jcp in Eq. (16) if Q = 0 at every point 
of sub-segment BC. 

In applying Eqs. (12) through (15) to a section MN 
of corrugated skin of constant thickness, the y’s can be 
measured to the median line of the corrugations, but 
Lyn must be the developed length of the corrugated 
skin instead of the length of its median line. 


CONCENTRATED AREAS 


It is customary to assume that the area of each longi- 
tudinal stiffener is concentrated at its center of gravity, 
the area of each stiffener on the compression side being 
augmented by the effective area of the skin associated 
with it. The work of shear stresses in such concen- 
trated areas is so small that it is usually neglected, the 
entire work of shear occurring in the skin of the outer 
perimeter and shear webs. These concentrated areas, 
however, do have an appreciable effect on the work of 
shear, because the transference to the skin of the effect 
of the statical moment of each concentrated area causes 
a change in the shear flow in the skin. Since the shear 
load in a stiffener is transmitted to the skin through its 
connecting rivets, the effect on the skin of the statical 
moment of each concentrated area, so far as shear flow 
in the skin is concerned, can be assumed to be concen- 
trated at a point on the skin coinciding with the center 
line of a single row of connecting rivets or midway be- 
tween the center lines of a double row of connecting 
rivets. 

When concentrated areas only are involved, it is con- 
venient to consider each concentrated area separately in 
the integral defining J (or 8), the total J (or 8) being 
the sum of these individual integrals. Thus, for the 
concentrated areas associated with a segment MN of 
skin, the J and 8 for bending about the x principal axis 


are: 


Jun = da; 9; Kyy (18) 
Bun = <a; y; Ajrw (19) 


in which a; is the area of any one of the concentrated 
areas associated with MN and 7’ is the point on MN to 
which the statical moment of a; is transferred. When 
the Q for a concentrated area is transferred to the skin 
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at the junction point MM of two adjacent segments 7 
and MN, it can be associated with either HM or MN 
but not with both of them. 


STRAIGHT-LINE SEGMENT OF CONSTANT THICKNESS 


Let MN be a straight-line segment of the skin of con- 
stant thickness fy;y, and let yy be the perpendicular 
distance from point P to line MN. If S is a point on 
MN and ds is the differential of arc along MN, dAys = 
dyn ds/2. Then, from the definitions of 8, J and A, 


Bun = Sit Ous bun ds/2 = dyn tun Jun/2 (20) 
Ayn = 6yn Lyy/2 = dyn tun Kun /2 (21) 


Note that Eqs. (2) and (21) are valid for concentrated 
areas as well as for skin effective in bending. 

For skin effective in bending, the substitution of any 
even integer for m in Eqs. (13) and (12) yield, re- 
spectively, 


Jun = (29a + Yn) L? yn /6 (22) 
Bun = (29), + Yn) Law tun Auwn/6 (23) 


The Jyy and Byy for the concentrated areas associ- 
ated with segment WN can be obtained from Eqs. (18) 
and (19). 


NEUTRAL AXIS 


If the skin is thin, the skin on the tension side of a 
cross-section is entirely effective in bending while on the 
compression side only the effective width of skin associ- 
ated with each longitudinal stiffener is effective in bend- 
ing. For this reason, the section properties of a cross- 
section cannot be determined until the location of the 
neutral axis (the line separating the tension and com- 
pression sides of the cross-section) is known. In 
general, the angle between the neutral axis and any 
orthogonal set of axes, u and v, through the center of 
gravity of the cross-section is a function of the section 
properties of the cross-section with respect to u and v 
and the components of the resultant moment about the 
u and v axes. Formulas for determining this angle are 
given in reference 3. Although the neutral axis of an 
airplane wing coincides with the beam principal axis 
(the one nearly perpendicular to the beam forces on the 
wing) only when the resultant shear force is actually 
parallel to the beam direction, the angle between the 
neutral axis and the beam principal axis is usually very 
small. Hence, it is generally accurate enough to as- 
sume that the neutral axis coincides with the beam 
principal axis when calculating the section properties of 
a wing cross-section. 


Axes Not PRINCIPAL ONES 


Let R, or R, stand for any one of the parentheses on 
the right side of Eqs. (5) through (11) according as there 
is bending about the x or y principal axes. Let x’ and 


y’ be any other pair of axes through the center of gravity 
of C,, let 6 be the angle (+ when counterclockwise) from 
the positive x’ axis to the positive x axis and let R,, and 
Ry be the same functions of y’ and x’, respectively, as 
R, is of y. Since R, contains the first power of y only, 
R,, R, and R, contain the first power only of x, y’ and 
x’, respectively. Moreover, since x = y’ sin @ + x’ cos @ 
and y = y’ cos @ — x’ sin @, 


Ry = Ry sin 6 + Ry cos 6 (24) 
R, = Ry cos @ — Ry sin 6 25) 
I, = Uy + Iy)/2 + Hyy/sin 26 (26) 
I, = de + Iy)/2 — Hyy/sin 20 (27) 
tan 20 = 2Hy»y,/(Ly — I,) (28) 


in which 7 stands for product of inertia. If one pre- 
fers, the values of R,, R,, J, and J, required in Eqs. (5) 
through (11) can be obtained directly from the x’ and y’ 
axes by means of Eqs. (24) through (27) without first 
having to find the coordinates of any point of the beam 
with respect to the x and y axes (the principal axes). 


EXAMPLE 


The cross-section shown in Fig. 3 will be taken as an 
example. The nose is a circle of 4-inch radius, but all 
other segments of the skin are straight lines. Segments 
ca and de are each 0.02 inch thick while segment dcbae 
is 0.03 inch thick. Only the portion of the skin below 
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Fic. 3. Two-cell wing used in the example. 


the neutral axis (the tension side) will be assumed to be 
fully effective in bending. The 7 upper concentrated 
areas, n, m, h, r, s, t and u each contain 0.12 sq.in., in- 
cluding the effective area of the associated skin. The 4 
lower concentrated areas, p, k, w and v, each contain 0.1 
sq.in. exclusive of any portion of the skin. It will be 
assumed that the Q of each concentrated area is trans- 
ferred to the skin as follows: areas h, k, u and v to webs 
ca and de at points coinciding with their respective cen- 
ters of gravity; areas m, m, r, s and ¢ to segment dc and 
dc extended at points directly opposite their respective 
centers of gravity; areas p and w to segment ea and ea 
extended at points directly opposite their respective 
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centers of gravity. By the usual methods, the center 
of gravity of Fig. 3 was located 0.36 inch below and 1.37 
inch to the right of the midpoint of line ca and the x 
principal axis (dj in Fig. 3) was rotated 8° 32’ counter- 
clockwise from line ae. It will be assumed that 07 is 
also the neutral axis of Fig. 3 and that point P is at the 
intersection of lines }j and ca. Table 1 gives the data 
for finding J and @ for the concentrated areas, the last 


TABLE 1 


Concentrated Areas 





Concentrated Areas 





Segment of Skin 


Name a y Name K A ayK ayA 
k 0.12 4.02 ha 375 0 180.90 0 
» 24 —2.90 ka 25 0 — 7.25 0 
ca 173.65 0 
m 0.12 4.32 m’'ba 618.88 52.56 320.83 27.25 
e @3 4.62 n'ba 552.21 47.99 306.15 26.61 
p 0.1 —2.30 p’a 133.33 6.86 — 30.67 — 1.58 
cha 596.31 52.28 
¢ 0.32 3.73 r'de 600 45.71 268.56 20.46 
S 0.12 3.43 s’de 1600/3 41.14 219.52 16.93 
t 0.32 3.14 t’de 1400/3 36.57 175.84 13.78 
u 0.12 2.84 ue 375 30.00 127.80 10.22 
v 0.1 —4.09 ve 25 2.00 — 10.22 — 0.82 
cde 781.50 60.57 
w 0.1 —3.50 we 400/3 6.86 — 46.67 — 2.40 
ae — 46.67 — 2.40 


two columns listing the quantities appearing in Eqs. 
(18) and (19). Table 2 gives data needed for Q, J and 
8 of segment ba of the skin when ba is divided into 6 


TABLE 2 


Segment ba of Skin 








Point y Cj (3-i/2)c; Aa A aiYi 
0 (b) 0 —1.26 —3.78 19.50 0 
1 —1.26 14.60 —18.40 
2 —2.22 —6.23 —12.46 10.50 —23.31 
3 —2.75 7.30 —20.08 
4 —2.97 —8.90 —8.90 4.87 —14.46 
5 —3.18 2.44 — 7.76 
6 (a) —3.39 —4.875 0 0 
Totals —21.265 —25.14 

Zsr — 260.50 





equal parts, the third, fourth and sixth columns listing 
quantities which appear in Eqs. (14), (13) and (12), 
respectively. Hence, for concentrated areas only (see 
Table 1): 


Jeq = 173.65, Jeng = 596.31, Ja = 781.50, 
Jee = —46.67, Besa = 52.28, Brae = 60.57, 
Bae = —2.40, Q, = 0.12(4.02 + 4.32 + 4.62) + 


0.1(—2.9 — 2.3) = 1.0352 


For skin alone (see Table 2 and Fig. 3): 
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Loa == 8.48, Lg — 3.43, Lie = 4.64, Lic = 
— 4.59, o. == 3.43, Die = § 


Ye = —3.39, y, = 

Jog = —25.14 X 8.48?/27 = —66.96 (Eq. (13)) 

Boo = —260.5 X 8.48 X 0.03/18 = —3.68 (Eq. (12)) 
Qra = —21.265 X 8.48 X 0.03/9 = —0.6011 (Eq. (14) 
Jp, = (2 X 0 — 3.39)3.43?/6 = —6.65 (Eq. (22)) 

Op, = 0.02 X 3.43 (—3.39)/2 = —0.1163 


Jj. = (2 X 0 — 4.59)4.642/6 = —16.47 (Eq. (22)) 
— 1.32 (Eq. (20)) 


8; = 8 X 0.02 (—16.47)/2 = 
Jue = {2(—3.39) — 4.59} 82/6 = —121.28 (Eq. (22) 
Bue = 3.43 X 0.03 (—121.28)/2 = —6.24 (Eq. (20)) 


For skin and concentrated areas: 


Jeo = 173.65 — 6.65 = 167.00 
Joa = 596.31 — 66.96 = 529.35 

Bsa = 52.28 — 3.68 = 48.60 

Jete = 781.50 — 16.47 = 765.03 

Bete = 60.57 — 1.32 = 59.25 

Joe = —46.67 — 121.28 = —167.95 

Be |= = 2.40 — 6.24 = —8.64 

Q, = 1.0352 — 0.6011 — 0.1163 = 0.3178 
QO:Kae = 0.3178 X 8/0.03 = 84.75 


Q,A,. = 0.3178 X 3.43 X 8/2 = 4.36 
Kc = 8/0.02 = 400, K,; = 8/0.02 + (8 + 4m)/ 
0.03 = 1085.55 
Kz = (8 + 8)/0.02 + (8 + 8)/0.03 = 1333.33 
A, =4X8+ 89 = 57.13,4. = 8X8 = 64 
I, = 22.66 
1085.55q, — 400g,’ = (529.35 — 167)/22.66 = 
15.99 (Eq. (9)) 
1333.33¢,, — 400g, = (84.75 + 167 — 
167.95 — 765.03)/22.66 = — 30.06 (Eq. (10)) 
Ge = 0.00722, g.’ = —0.02038 
Dp/2 = 57.13(0.00722) + 64(—0.02038) + 
(59.25 — 48.6 + 8.64 — 4.36)/22.66 
Dp = —0.47 inch (to left of P) 


(Eq. (11)) 


The Dp for bending about the y axis is found by a 
similar process with J, replacing J, and x replacing y in 
the formulas for Q, J and 8. An actual airplane wing 
would not be much more difficult than this example, the 
only added complication being that ae would be a curved 
line and would require an additional table similar to 
Table 2 for the calculation of its J and B. 
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SUMMARY 


It is shown how the thrust and the thermal efficiency of a 
ducted radiator vary with the design of the duct, the amount of 
radiator heating and the airplane speed. The type of duct 
commonly used for drag reduction is found to serve well for thrust 
recovery. Nevertheless, the thermal efficiency of the arrange- 
ment is so low that the thrust obtainable is unlikely to offset more 
than about one-third of the drag of the radiator and duct. 

A study is made of the thrust given by a jet of fluid released ina 
duct. Conditions are found to exist under which the thrust of 
the ducted jet is greater than that of the jet unducted. If the 
jet consists of cool fluid, the thrust augmentation is secured by a 
venturi-shaped duct. With a hot jet, an expanding-contracting 
duct is required, as for a propulsive radiator. The analysis sug- 
gests that the potential benefits of exhaust ducting are of the 
same order as those of radiator ducting. 


INTRODUCTION 


HE power delivered by a moving engine can be in- 
creased through a utilization of the cooling air and 
the exhaust gases as propulsive jets. 

That a propulsive effect should be obtainable from 
the cooling air was suggested several years ago by 
Meredith.' Since then, various writers? * have elabo- 
rated on the idea and clarified the basis of its working. 
The cardinal point is that the expansion of a stream 
of fluid under heating will result in a speeding up of the 
fluid, and hence in a propulsive reaction, only if the 
stream is confined between properly shaped boundaries. 
The variation of the propulsive effect with the charac- 
ter of these boundaries is examined from an elementary 
theoretical viewpoint in Part I of this paper. 

The possibility of deriving a thrust from the exhaust 
gases by merely pointing the exhaust stacks backward 
has long been recognized. The propulsive force that 
can be secured in this way is, however, relatively small.® 
A greater force could, in principle, be produced through 
a diffusion of the kinetic energy of the exhaust jets 
into a larger mass of fluid. This has led to the proposal 
of thrust augmentors in the form of various types of 
ducts around the exhaust stacks.7"" An elementary 
theory of the performance of such devices is developed 
in Part II of the paper. 


Part I 
RADIATOR PROPULSION 


Radiators in Ducts 
The cooling of aircraft engines requires a dissipation 
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of heat into the atmospheric air. The bodies through 
which this dissipation takes place are termed radiators. 
If the engine is air-cooled, the engine itself is the radia- 
tor. When the engine is cooled by liquid, most of the 
heat is disposed of through a radiator set apart from 
the engine. In either case, the radiator is customarily 
enclosed in a cowling or duct that provides an expand- 
ing air channel in front of the radiator and a contract- 
ing channel behind the radiator. The primary pur- 
pose of these fairings is to keep down the drag by guid- 
ing the air through and past the radiator with a mini- 
mum of disturbance. As it happens, however, the 
same type of fairing is capable of extracting propulsive 
energy from the cooling air. The net drag reduction 
achieved by a fairing thus varies with the extent to 
which this propulsive effect is utilized. 


Thermodynamics of Propulsive Action of Ducted Radiators 


The derivation of a propulsive effect from the heating 
of a ducted stream can be explained in terms of simple 
thermodynamics if the fluid is a gas, like the air in the 
problem under study: 

While slowing down and gaining pressure in its mo- 
tion through the entrance section of the duct, the air 
compresses essentially adiabatically. Then heat is 
added at approximately constant pressure as the air 
passes through the radiator core. Next, the air under- 
goes essentially adiabatic expansion as it speeds up and 
loses pressure in the exit section of the duct; and finally 
it is cooled at atmospheric pressure until it returns to 
the condition of the entering air. Through this ther- 
modynamic cycle, heat energy is changed into me- 
chanical energy with an efficiency 7 = 1 — (1/r)*~}, 
if r signifies the ‘compression ratio”’ of the cycle (7.e., 
the ratio of the density of the air at the radiator to the 
density of atmospheric air), and & is the ratio between 
the specific heats of the cooling air at constant pres- 
sure and at constant volume. 

If the process is viewed from a system fixed in the 
radiator, so that the air moves while its constraints are 
stationary, any work absorbed by the air must ulti- 
mately appear as an addition to the air’s kinetic energy, 
and so result in an increase in the momentum of the 
ducted stream. The creation of these energy and 
momentum changes requires a backward push on the 
air by the radiator and the duct. The air, in turn, 
reacts on the radiator and duct with an equal propulsive 
force. 
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To determine the useful work accomplished, it is 
necessary to examine this process from a system fixed in 
the atmosphere. The radiator then advances at the 
velocity of the airplane, while the cooling air stirs 
relatively little. Most of the work done by the mutual 
forces between radiator and air thus is done on the 
radiator, and is therefore utilized for propulsion. The 
work not utilized is found in the final kinetic energy of 
the cooling stream relative to still air. To arrive at 
the propulsive efficiency of the process, it is thus neces- 
sary to multiply the thermal efficiency by a factor 
analogous to the Froude efficiency of a propeller. 

Because of the uncertain value of this Froude factor, 
and the inaccuracy of the assumption that the heating 
of the air takes place at constant pressure, the foregoing 
analysis does not yield precise results. A number of 
general conclusions can, however, be drawn from it. 
For example, the efficiency of the radiator as a pro- 
pelling device is seen to increase rapidly with the air- 
plane speed, because the compression of the duct air, 
which determines the thermal part of the efficiency, 
varies as the square of the speed. The compression 
depends also on the ratio of the area of the radiator core 
to the entrance area of the duct. The maximum effect 
would be secured if this ratio were infinite, so that the 
flow would reach its full stagnation pressure at the core. 
An effect almost as great results, however, if the core 
area is but three or four times the entrance area, as in 
ducts actually used in practice. A parallel-sided duct 
would have no compressive effect and would therefore 
have zero efficiency. This is in line with the observa- 
tion that no pressure forces could be exerted on the 
walls of such a duct in the direction of the motion. A 
venturi-shaped duct would have a negative efficiency, 
and so would produce a drag under heating. 

A rough estimate of the magnitude of the propulsive 
effect of a well-ducted radiator may be made by as- 
suming that the full stagnation pressure is realized at 
the radiator core, and that the kinetic energy losses are 
negligible. Excessive as such an estimate is, it indi- 
cates propulsive efficiencies of only a few per cent at 
the highest airplane velocities yet reached. The 
power obtainable by radiator ducting thus is unlikely 
to be more than a small fraction of the propulsive power 
of the engine itself. 


Hydrodynamic View of Propulsive Effect of Radiator 

Ducting 

The thermodynamic process discussed in the pre- 
ceding Section is but one phase of the energy trans- 
formations taking place in the ducted flow. In a 
broader view of the problem, thermodynamics merges 
into fluid dynamics and is covered by the latter. Be- 
cause of the smallness of the density variations in the 
flow at ordinary speeds, the flow is esserttially as it 
would be in the total absence of compressibility effects. 
With allowance for a density change brought about by 
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the heating, the problem can thus be treated as one in 
hydrodynamics. 

The cross-sectional area A, velocity v and pressure p 
of the cooling stream at stations far ahead of the radia- 
tor, immediately in front of the core, immediately in 
back of the core, and far behind the radiator will be 
marked by the indices shown in Fig. 1. With p; 
the density of the unheated fluid, and 2 the density of 
the heated fluid, Bernoulli’s equation applied in front 
the radiator and behind the radiator furnishes the 


relations 


1/sp:017 + pi = '/opite? + pe (1) 
and 
1/sp203" + ps = '/apevs? + pa (2) 
Continuity requirements make 
piAyw: = pidote = prAos = prAsrs (3) 
The momentum equation for the flow through the core 
is 
(b2 — ps)Az = prteA2(vs — v2) (4) 
and for the thrust T 
T = prAove(vy — 2) (5) 


From these seven equations, it is possible to deter- 
mine seven of the variables in terms of the remaining 
variables. If the former are taken as po, v2, 3, v3, As, % 
and 7, the solution for T proceeds through the equations 








‘a = pains ad 1) (6) 
proAs 
and 
a-Vb+@]-@y ° 
— = i 1 a OE pe (7) 
ie ak 2) A: 
to 








T = pAw,? of =) - (24) - ; . 
EOE 


Of the variables entering into this expression, p, and 
v are given directly. 2 is determined by the tempera- 
ture rise allowed in the cooling stream. This tempera- 
ture rise, in conjunction with the specified rate of heat 
dissipation, fixes the mass flow p:A,v, through the radia- 
tor, and so determines A;. From the mass flow and 
the heat dissipation, it is finally possible to determine 
A; for any given type of radiator core. 

The quantity 7/p,A.v,*, representing a thrust co- 
efficient based on the initial area and velocity of the 
cooling stream, is plotted in Fig. 1. The curves show 
that a positive thrust is obtained only when the ratio 
A,/Az lies below a critical value which is always less 
than unity and decreases as the heating increases. 
For any given amount of heating, the thrust is greatest 
when A,/A2 equals zero, 7.e., when the radiator core 
area becomes infinite. This result agrees with the con- 
clusion reached in the thermodynamic analysis before. 
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Thrust Calculation Based on Exit Area of Duct 


If the duct extended far enough forward and back- 
ward, and approached constant cross-sections toward 
both ends, the cooling fluid would enter the duct at the 
undisturbed velocity v;, and the cross-sectional area of 
the mouth of the duct would be A;—assuming that the 
areas of mouth and exit were in the ratio A,/A; called 
for by Eq. (7). The various quantities entering into 
the thrust equation could then be computed directly 
from the layout of the duct. The problem is, however, 
complicated by the fact that the proper ratio A,/A; 
varies with the ratio p2/p:, and hence with the heating. 
Because of this fact, the entrance and exit openings 
of a fixed duct can be in the correct ratio only under 
one heating condition. When the ratio between the 
duct openings is incorrect, the proper ratio A,/A3; 
must establish itself through changes in the cross- 
section of the cooling stream beyond the ends of the 
duct. The identification of the entrance opening with 
A, and of the exit opening with A; is then no longer 
permissible. 

For a rigorous determination of the actual areas A; 
and A; under these circumstances, it would be neces- 
sary to resort to experiment. An approximate ra- 
tional method of obtaining the areas can, however, be 
based on the known fact that the fluid in a jet emanating 
from a parallel-sided pipe is roughly at the pressure of 
the fluid surrounding the jet. This fact suggests that 
the flow through a duct with a parallel-sided end sec- 
tion will always be such as to yield a discharge jet 
of nearly constant cross-sectional area. A; may, 
accordingly, be taken to equal the exit area of the duct, 
while A; is left to assume whatever value it must have 
to satisfy Eq. (7). 

Through this argument, the details of the entrance 
section become irrelevant—so long as the section is wide 
enough to permit an adequate inflow. Todo away with 
the restriction to a parallel-ended exit section, it is 
necessary to discover how the ultimate cross-sectional 
area of a contracting or expanding jet is related to the 
jet area at a point of given convergence or divergence. 
Some light on this question can be obtained from 
studies of propeller slipstreams.’ The indications are 
that ducts with walls converging or diverging toward 
the exit at angles up to 10 or 15 degrees will yield jets 
whose ultimate cross-sectional area differs from the 
exit area by but a few percent. A; may thus be taken 
to equal the exit area of the duct when the duct end 
is tapered, as well as when it has parallel sides. 

The area A), in terms of which the thrust was ex- 
pressed by Eq. (8), is related to A; through Eq. (7). 
Solved for the ratio A,/A_ instead of A/A3;, the latter 
equation gives 


~ 42 1+ (4) | -1 ) 
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Fic. 1. Dependence of radiator thrust coefficient 
T /p,A,v;? on ratios p;/p, and A;/Ag. 


The thrust equation can therefore be written explicitly 
in terms of A; as 


Pl 
p2 


vile me A;\* A;\* 7 
(e|1 (3) | . (3) 
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For the ratio p;/p2, it is possible to write 











m= 1+7(As) (11) 
p2 
with 7 the coefficient of thermal expansion of the cooling 
air at essentially constant pressure, and A? the differ- 
ence between the temperatures of the heated and un- 
heated air. This temperature difference is deter- 
mined by the equation 


H 

ce. ae 32.2 pA, 73c( At) (12) 

778 
in which H is the rate of dissipation of heat energy, in 
foot-pounds per second, and ¢ is the specific heat of the 
cooling fluid, in B.t.u. per pound per degree Fahrenheit. 
The equation for p;/p2 can thus be written 





ioe, 1s Hr (13 
ios (32.2) (778) pA wie 


or, with 7/(32.2)(778) c lumped into a constant K, 


_H_ (As As 


KpiAwi { As ; Aj 





m1 KpiAwy a 


H 2 
ae A;\? A;\? (14) 
i+ KpAw ye [3 * (4) | - (2) 


or, after transposition, squaring and reduction, 


a 
pa 
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H H A;3\" 
pl = a = al =o a 
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This expression, through Eq. (10), now yields the thrust 
directly in terms of the independent variables p, 1, A2, 
A;, H and K. 

The dissipation of the heat energy H enables the radi- 
ator and duct to do the useful work 77. The arrange- 
ment may therefore be said to have an efficiency 7 = 
Tv,/H. Thus, with the values of 7 and H given by 
Eqs. (10) and (14), 


pi / LA N\2 1] /A.\2 
(a NalL+ Ga) I- Ge 
_ <a p2 : As eer (16) 
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For given ratios p;/p: and A;/Az, the efficiency is seen 
to be proportional to the square of the velocity 1. 
This result is in agreement with the conclusions reached 
through the thermodynamic analysis before. The 
plot of Eq. (16) in Fig. 2 also confirms the expectation 


that the efficiency would be greatest with an infinite 
core area, but that the efficiency is only slightly less 
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Variation of thermal efficiency of ducted radiator 
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with core areas equal to three or four times the duct 
openings. 

The heat power dissipated by the radiator of a gaso- 
line engine normally equals about 50 per cent of the 
engine’s indicated power, while the propeller delivers 
in useful form about 75 per cent of this indicated power. 
The ratio of the energy dissipated through the radiator 
to the useful work done by the propeller is therefore 
about */;. Multiplication of this ratio by the thermal 
efficiency of the radiator gives the ratio of the propulsive 
power of the radiator to the thrust power of the propel- 
ler, or, what is the same thing, the ratio of the radiator 
thrust 7 to the propeller thrust 7,. Fig. 3, based on 
the assumption that c = '/, and 7 = '/529, so that K = 
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3,250,000, shows the relation between 7 and 7, for a 
typical case (p:1/p2 = 1.2, A2/A3; = 4). 


Effects of Friction Losses in Radiator and Duct 


The assumption of an ideal flow in the foregoing 
analysis helps to bring out the working essentials of 
radiator propulsion, but leaves the solution of the prob- 
lem incomplete, especially through the neglect of the 
friction loss across the radiator core. To allow for this 
loss, it is necessary to add a quantity A» App to the right- 
hand side of Eq. (4). The thrust equation (8) then 
becomes 


T = pAgt y" ay a | e (mar 
prAwy { ar + (4: re coe 
_ i} (17) 


Because the ratios p:/p2, A1/A2 and Apo/!/2p,2;? all 
influence each other in a manner peculiar to each case, a 
quantitative discussion of this modified thrust equation 
will not be attempted here. Qualitatively, it is evident 
that the presence of the term Apo/!/2p11;? reduces the 
radical in the expression, and hence cuts down the 
thrust. As this radical represents the quantity 
p1A1/p2A3, it appears also that the exit area of the duct 
must be increased if the values of p;, p2 and A; are to be 
the same as before, 7.e., if the cooling capacity of an 
actual system is to come up to that of a frictionless one. 
The whole analysis presumes, of course, that the effec- 
tive core area is not reduced through flow separation 
in the entrance section. 

Besides the friction loss in the radiator itself, there 
is also a loss from the air resistance of the duct. Ex- 
perimental indications are that the combined losses 
always overbalance the gain to be had from the heat- 
ing. For density ratios up to 1.2, for example, as 
found in conventional installations, Fig. 1 shows a 
thrust coefficient of less than 0.10, while the drag co- 
efficient of a typical unheated radiator and duct sys- 
tem is about 0.30 when based on the same area and 
By raising the density ratio to a higher 
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figure, it would be possible to secure a higher thrust 
coefficient. This gain would, however, be largely offset 
by the increased friction losses resulting from the greater 
viscosity of the hotter air and from the extra cooling 
surface needed to make up for the reduced tempera- 
ture difference between the radiator and the air. A 
thrust somewhat above one-third of the gross drag of 
the warm radiator and duct thus seems the best that 
may be expected on the basis of present knowledge. 


Part II 
EXHAUST PROPULSION 
Ducted Jets 


The kinetic energy mv*/2 of a mass m with a velocity 
v is connected with the momentum mv through the 
identities 


mv?/2 = (mv)(v/2) = (mv)?/2m 


The momentum corresponding to a given kinetic energy 
thus is relatively great when the energy is distributed 
through a large mass moving slowly, and the momen- 
tum is relatively small when the energy is concentrated 
in a little mass moving fast. It is then possible, in 
theory, to secure an increase of momentum without 
expenditure of kinetic energy through a dissemination 
of the energy into a larger mass. 

This theoretical possibility proves difficult to realize 
in practice. Merely mixing a high speed jet with 
slower moving fluid in an open stream, for example, 
yields no benefit, because so much energy is lost in the 
transfer from jet to fluid that the increased mass with 
the residual energy has only just the same momentum 
as the original jet. That no momentum change could 
be produced by this kind of mixing follows, of course, 
directly from the absence of any external agent capable 
of furnishing the force necessary for such a change. 
The problem of thrust augmentation thus is the dual 
one of reducing the mixing losses and of providing a 
support for an accelerating force. This dual require- 
ment can be satisfied by a non-parallel-sided duct 
around the jet. 


Hydrodynamic Theory of Ducted Jets 


In the ordinary exhaust collector, which also serves 
as a muffler, the exhaust gases are slowed down well be- 
low the velocity of sound. Compressibility effects are 
therefore of secondary magnitude, and the flow may be 
studied by the hydrodynamic method already used in 
connection with the radiator problem. 

The setup considered is shown in Fig. 4. Exhaust gas 
with a density p, and velocity v, is assumed to be dis- 
charged into a duct from a parallel-walled pipe of cross- 
sectional area A,. In accordance with the principle 
noted in the study of radiator jets, the gas is assumed 
to emerge from the pipe at the pressure /2 existing in 
the duct around the mouth of the pipe. The outer 
stream has a density p; and a velocity 7 far ahead of 


the duct. The portion of the stream passing through 
the duct has a cross-sectional A; far ahead, and a net 
area A» about the mouth of the inner pipe, past which 
this ducted flow has a velocity v2. Mixture of the ex- 
haust gas and of the augmentor fluid finally yields 
a jet of density p3, cross-sectional area A, and velocity 
vs at the open stream pressure /; far downstream from 
the duct. 

The mixing of two flows in a channel of variable cross- 
section is a process that cannot be treated by elementary 
methods. In the present analysis, it is therefore as- 
sumed that the mixing takes place in a parallel-sided 
section of the duct. Since the duct walls then can 
exert no axial force on the fluid in this mixing section, 
the momentum changes occurring in the section are 
determined by the pressures, densities and velocities 
at the ends of the section. The effects of the mixing 
can, accordingly, be allowed for in the same general 
manner as the effects of the heating in a radiator core. 

The results obtained through the assumption of a 
cylindrical mixing section hold equally when the mixing 
takes place in an infinitesimal length of a duct of vari- 
able cross-section. Such instantaneous mixing could 
be approximated in practice by the use of a large num- 
ber of jets instead of one central jet. When the mixing 
process extends over a distance in which the duct sec- 
tion varies perceptibly, as it usually does under actual 
conditions, the theory does not apply with rigor. It 
should, however, still remain a fair approximation so 
long as the gases are uniformly mixed by the time they 
finally leave the duct. 

Bernoulli’s equation applied in front of the mixing 
section and behind the mixing section gives the rela- 
tions 


1/spi01? + pi = */opiv2* + po (18) 
and 
1/spq3" + ps = '/opsve? + pr (19) 
while continuity requires that 
PAW, = prAgr (20) 
and 
piAgwe + PeA.v, = ps(A2 + A,)vs3 = p3Aq; (21) 


if v3 is the uniform velocity of the mixed stream at the 
end of the mixing section. The momentum equation 
for the flow through the mixing section is 


(p2 — ps)(A2 + As) = p3(A2 + A,)v5? — 
piAgve” — p,A,v,* (22) 
if p; denotes the pressure of the mixed fluid at the ve- 
locity v3. Finally, there is the thrust equation 
T = pAge — pA? (23) 


Through these seven equations, it is possible to express 
seven of the variables in terms of the others. If the 
dependent variables are taken as po, v2, ps, Us, As, V4 
and 7, the solution leads through the relations 
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.. aha, + p1A qv (24) temperature, a duct with a throat area about 0.7 times 
i‘ p3A4 iS the inlet area would give the maximum improvement, 
and but even so would only raise the thrust about 5 per 
A os pA ,v, + prA 
. =  —— —————————————————— oe — — — —<— = (or 
nye|s — (4) ]— al etactnden y Bedod Ander 
Ps 2 ps A;+A, px(A2+A,)  p3A2(A2 + A,) 
to 
_ 9 Pp, Av 
T = piAyq,? ! — + | 
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The complexity of this expression makes it difficult to 
recognize the influences of the different basic variables 
on the thrust. Some idea of those influences may, how- 
ever, be gained by assigning reasonable values to some 
of the quantities, and finding how the thrust varies 
with the remaining one or two quantities. For ex- 
ample, a jet of constant mass flow and speed may be 
assumed. This fixes the thrust that would be given 
by the jet if released in free air. Assuming, next, that 
the exhaust gas is to mix with air in a constant propor- 
tion, the entrance area A, the basic stream velocity 
v,; and the basic stream density p, may be assigned 
fixed values. This leaves to be varied the density p, 
of the exhaust gas, and the area A; + A, of the parallel- 
sided portion of the duct. 

Fig. 4 shows the result of this procedure for A, = 
1 sq.ft., 1 = 200 ft./sec., p, = 0.00239 slugs/cu.ft., 
and an exhaust mass flow p,A,v, = 0.03 slugs/sec., 
which corresponds to the air capacity of a 500 hp. en- 
gine. The exhaust velocity is taken as 600 ft./sec., 
about what would be used in practice. The exhaust 
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Fic. 4. Effect of ducting on thrust of exhaust jet. 
flow would then require a pipe 1.96 inches in diameter 
if the flow consisted of cold air, or 4.2 inches for exhaust 
gas at 2000°F.—the properties of the exhaust gas being 
taken as those of air, since both the gas and the air 
contain mostly nitrogen. : 

Under the conditions assumed, the exhaust thrust in 
free air would be 18 Ibs. If the exhaust were at room 
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cent. As the throat area of the duct is made less than 
the entrance area, the thrust becomes less than that 
of the free jet. To secure a thrust augmentor effect, 
it is thus necessary to speed up the duct air so that it 
will mix with the exhaust without too much kinetic 
loss. If this speeding up is carried too far, however, 
the net throat area becomes so small that a very 
large exit area is necessary to suck through the desired 
amount of duct air. This slows down the emerging 
stream, and results in a diminished, or even negative, 
thrust. 

From the other curves in Fig. 4, it is seen that mild 
warming of the exhaust reduces the attainable thrust 
augmentation, but that an increased augmentation 
again is possible with a very hot exhaust. The aug- 
mentation in the latter case requires a duct with a 
throat area larger than the inlet area. This means 
that the augmentor action has become predominantly 
the result of a utilization of the exhaust heat, analogous 
to the utilization of the heat furnished by a ducted 
radiator. Asin the case of the radiator, a large amount 
of heat is needed for any appreciable effect. Thus, a 
maximum thrust of 22 lbs., or an augmentation of but 
4 Ibs., could be obtained if the temperature of the 
exhaust gas were 2000°F., 7.e., about 300°F. above the 
normal temperature of such gases. 

These observations suggest that it is impossible to 
utilize for thrust augmentation both the kinetic and 
the thermal energy of the exhaust. In a venturi- 
shaped duct, which serves to extract kinetic energy, 
there is a loss for thermodynamic reasons. A duct of 
the expanding-contracting type gives beneficial re- 
sults on the thermodynamic side, but produces in- 
creased kinetic losses; and though the extra heating 
brought about by the losses in the latter case may be 
indirectly helpful, it has very little effect because of 
the low efficiency of the thermodynamic process. Any 
marked thrust augmentation presupposes, then, a large 
disparity between the thermal and kinetic energies 
extractable from the exhaust, with a consequent possi- 
bility of drawing on one of these energies without serious 


losses on account of the other. 
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The full possibilities of thrust augmentation are, of 
course, not indicated by Fig. 4, which merely represents 
one special case. Better results could certainly be 
found at higher velocities ,, where the propulsive ef- 
fect of the heating would be more pronounced. As the 
heat energy in the exhaust jet is about four times as 
great as the heat energy passing through the radiator, 
the exhaust thrust might be raised to something like 
four times the radiator thrust by proper ducting. 
This would represent a substantial augmentation under 
favorable conditions. If the radiator is already ducted, 
a simple way of utilizing the exhaust would then be 
to discharge it in the radiator duct just behind the radia- 
tor. In an air-cooled installation, this would be ac- 
complished by releasing the exhaust inside the cowling 
without other manifolding. 


Results Based on Exit Area of Duct 


The area A;, on which the thrust has been based in 
the analysis so far, must itself be determined through a 
supplementary calculation. To this end, it is con- 
venient to put Eq. (25) into the form 


[mn fl: 27] +of88 + adap} 
+2[ nam] ode. [424+ oA" ]} + 
[nen [48+ oa] -[2T- 
; 
oe Proeyt =0 (27) 


and to introduce the abbreviations 
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through which the equation is reduced to 


(u — A)(p1Ai01)? + 24 (p1A101) (pA 2) + (u — v) 
(p,A,v,)? _ 0 (28) 
Then, 
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and, by Eqs. (23) and (24), 
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with 7, = p,A,v,? the thrust that would be furnished 


by the exhaust jet if there were no duct. The quantity 
in curved brackets might thus be termed the augmenta- 
tion factor for the ducted arrangement. 

Instead of the area A; used before, the thrust expres- 
sion now contains the area Ay. By the logic already 
followed in the case of the ducted radiator, A, may be 
assumed to equal the exit area of the duct if the expan- 
sion or contraction of the duct toward that area is rea- 
sonably gradual. This leaves the thrust completely 
defined by quantities known directly in any given case. 

Aside from making the solution of the problem more 
definite, the use of A, instead of A; has the advantage 
of removing the difficulty of handling the case of v, = 0, 
for which the thrust expression in terms of A;, with 
A, = ©, would become indeterminate, while the 
formula based on A, gives at once 
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Experimental Sidelights on Augmentor Theory 


In a few rough preliminary tests on ducted jets, the 
results have been found to agree closely with those 
predicted by the foregoing theory.® Earlier findings 
by Schubauer”™ suggest that the cool jet augmentation 
shown in Fig. 4 for a particular case is typical of the 
thrust increases that may be secured without heating. 
The more pronounced augmentor effects found by 
Jacobs and Shoemaker"! are also in line with the theory. 
For, since the driving jet in the experiments of these 
authors came from a high pressure tank at nearly at- 
mospheric temperature, the jet was presumably much 
colder than the augmentor air. In the case of a duct 
of the venturi type, such as was used, cooling produces 
the same sort of propulsive effects as heating does in 
the case of an expanding-contracting duct. The direct 
augmentor ability of the duct was therefore supple- 
mented by helpful thermodynamic effects, and a sub- 
stantial total augmentation seems not unreasonable. 


Exhaust Collectors vs. Individual Stacks 


The foregoing analysis presumes steady flow condi- 
tions and applies therefore only to exhaust streams fed 
by numerous cylinders. Such collected streams are 
generally much slower than the free blasts from cylin- 
ders with individual stacks. This means that exhaust 
collection reduces the net impulse /Fdt on the dis- 
charged gases, and hence the work {Fds = /{Fo,dit = 
v, { Fdt done by the net exhaust reaction F on the en- 
gine moving at the velocity v,. If the noise and glare 
of unmuffled exhaust flashes could be tolerated, or 
could be blanketed from view and hearing, it would 
thus be advantageous to use individual stacks, prefer- 
ably shaped into expanding nozzles. There is some 
indication® that augmentor effects can be secured with 
intermittent jets as well as with continuous ones. It 
seems, however, hopeless to attempt a treatment of 
these effects along the elementary lines of this paper. 
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HONORARY FELLOW AND FELLOWS ELECTED 


The Fellows of the Institute have elected Dr. George W. Lewis 
an Honorary Fellow. Dr. Lewis is Director of Research of the 
National Advisory Committee for Aeronautics and a Past 
President of the Institute. He received the Daniel Guggenheim 
Medal in 1936 and delivered the Wilbur Wright Memorial Lec- 
ture before the Royal Aeronautical Society of Great Britain in 
1939. 

Ten American members of the Institute were elected as pro- 
vided in the By-Laws to the grade of Fellow in 1940. The new 
Fellows are: 

Allan Chilton, Project Engineer, Wright Aeronautical Corp. 
Charles H. Colvin, Acting Chief, Instruments Div., U.S. Weather 

Bureau. 

Smith J. De France, Engineer in Charge, Ames Aeronautical 

Laboratory, N.A.C.A. 

Melvin N. Gough, Senior Airplane Test Pilot, Langley Memorial 

Aeronautical Laboratory, N.A.C.A. 

Sam D. Heron, Research Engineer, Ethyl Gasoline Corp. 
Paul H. Kemmer, Materiel Div., U.S. Army Air Corps, Wright 

Field. 

Paul Kollsman, Director & Vice-Pres., Square D Company, in 

charge Kollsman Instrument Div. 

George A. Page, Jr., Chief Engineer, St. Louis Airplane Div., 

Curtiss-Wright Corp. 

F. W. Reichelderfer, Chief, U.S. Weather Bureau. 
Igor Alexis Sikorsky, Aero. Engineer, Vought-Sikorsky Aircraft, 

United Aircraft Corp. 


Dr. Louis H. BAUER RECEIVES JOHN JEFFRIES AWARD 
FOR 1940 


Louis H. Bauer, M.D., of Hempstead, N.Y., was named the 
first recipient of the John Jeffries Award recently established by 
the Institute, for his pioneering work and continuing activity in 
advancing the interests of aviation medicine through twenty-one 
years of teaching, research, organization, editing and contribu- 
tions to the literature. In recognition of the importance to avia- 
tion of scientific endeavor in the field of medicine, the Institute 
presents the John Jeffries Award annually for outstanding con- 
tributions to the advancement of aeronautics through medical 
research. 

Dr. Bauer was graduated from Harvard College and Medical 
School and entered the Medical Corps of the U.S. Army in 1913. 
He was the first commandant of the School of Aviation Medicine 
of the Army Air Corps which he organized at Mitchel Field in 


1919. He was the first Medical Director of the Aeronautics 
Branck of the Department of Commerce. Dr. Bauer organized 
and was the first president of the Aero Medical Association and is 
now editor-in-chief of its official publication, the Journal of Avia- 
tion Medicine. He is chief cardiologist at Meadowbrook Hospi- 
tal, Hempstead, N.Y., a member of the staffs of other Long Island 
hospitals, and consultant in aviation medicine and cardiology to 
the Civil Aeronautics Administration. He is a Fellow of the 
American College of Physicians. 


HENRY G. HOUGHTON NAMED FOR ROBERT M. LOSEY 
AWARD FoR 1940 


Henry G. Houghton, Jr., Assistant Professor of Meteorology at 
Massachusetts Institute of Technology, was chosen by the Com- 
mittee of Award to receive the 1940 Robert M. Losey Award for 
outstanding contributions to the science of meteorology as ap- 
plied to aeronautics, presented for the first time at the Honors 
Night Dinner on January 28th. He receives the award for 
fundamental research in physical meteorology on the processes of 
condensation in the atmosphere. 

Mr. Houghton was born in New York City in 1905. He brings 
to the field of meteorology the training and research methods of a 
physicist and electrical engineer, having studied at Drexel Insti- 
tute of Technology and received the degree of S.M. in Electrical 
Engineering from Massachusetts Institute of Technology in 1927. 
Following graduation, his meteorological work at the Round Hill 
Research Station and Massachusetts Institute of Technology 
dealt with the fundamental physical properties of atmospheric 
fog, apparatus and methods for measuring fog particles and the 
resistance of fog and clouds to the transmission of visible light 
and radiant rays. This led to his development, with the assist- 
ance of W. H. Radford in 1938, of a means for clearing up fog 
over local areas by the use of a calcium chloride spray. This 
method, described by Dr. Sverre Petterssen in his recent Wright 
Brothers Lecture, is one of the first practically tested ways for 
artificially dissipating fog, although its application is limited. 

Mr. Houghton’s recent research has been an extension of simi- 
lar studies in atmospheric condensation processes; i.e., the forma- 
tion of fog, clouds, water droplets and ice particles in clouds and 
icing on aircraft in flight. His work on the underlying principles 
of these phenomena and precise measurement of their characteris- 
tics have contributed substantially to greater exactness in the 
science of weather forecasting and increased safety in air travel. 
He is a Fellow of the American Association for the Advancement 
of Science and a member of the American Physical Society, the 
American Meteorological Society, Sigma Xi and Tau Beta Pi. 
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REPORT ON NINTH ANNUAL MEETING TO BE IN MARCH 
ISSUE OF JOURNAL 


Since this issue of the Journal went to press before the Honors 
Night Dinner and the Ninth Annual Meeting of the Institute, 
January 28th to 3lst, complete accounts of the Dinner and Meet- 
ing will be published in the March issue. 


PERSONNEL OPPORTUNITIES 


The Personnel Bureau serves individual members, as well as 
organizations seeking to employ aeronautical specialists. Any 
member or organization may have requirements listed without 
charge. 


Wanted 

Small California aircraft manufacturing company with a ply- 
wood plastic project needs several senior layout men with stress 
analysis experience; also a project engineer. Please include com- 
plete information regarding experience and salary. Address 
reply to Box 121, Institute of the Aeronautical Sciences. 

Aircraft Accessories Sales Engineers. Positions open to young 
men of proven ability in this field or allied fields. No applications 
considered without complete historical background and proof of 
ability. Must be willing to travel and locate anywhere in U.S. 
Address reply to Box 122, Institute of the Aeronautical Sciences. 

Aircraft Accessories Design or Project Engineers. Positions 
open to young men of proven ability in this field or allied fields. 
No applications considered without complete historical back- 


ground and proof of ability. Location Cleveland. Address 
reply to Box 123, Institute of the Aeronautical Sciences. 


Available 
Executive, associated with aircraft industry for 24 years. 
B.S.E. (Ae.E.) and Ae.E. degrees. Aero. Engineer, U.S. Army 
and U.S. Navy. Consulting Engineer, Project Engineer and 
Chief Engineer. Address reply to Box 124, Institute of the Aero- 
nautical Sciences. 


STUDENT BRANCHES 


University of California. Mr. A. F. Bonnalie, of United Air 
Lines Transport Corporation, gave a talk on ‘“‘Economic High- 
lights in Air Transportation” at a meeting held on October 3rd. 
A motion picture describing the Langley Memorial Aeronautical 
Laboratory of the N.A.C.A. was shown and discussed by Mr. 
Carlton Bioletti at a meeting held on November 7th. 

Carnegie Institute of Technology. The Fourth Wright 
Brothers Lecture by Dr. Sverre Petterssen on ‘‘Recent Fog In- 
vestigations” was delivered by Walter E. Mortlock at a meeting 
held on December 19th. D. M. O’Keefe, Chief Meteorologist, 
of Pennsylvania-Central Airlines, and his assistant, O. P. Braun, 
participated in a discussion which followed. 

Oregon State College. At a meeting held on December 11th, 
R. E. Herr, of the General Inspection Division of the C.A.A., 
spoke on the engineering aspects of his work. 

Rensselaer Polytechnic Institute. On December 5th, Dr. Paul 
Hemke, Head of the Aeronautical Engineering Department, gave 
a talk on “High Speed Visual Flow Research.” 
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of the AERONAUTICAL SCIENCES 


The Institute of the Aeronautical Sciences invites 
both members and non-members from any country to 
submit papers for publication in the Journal of the 
Aeronautical Sciences. The Journal, following the 
practice of other scientific publications, does not pay for 
contributions. 


The following directions for the preparation of papers, 
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more than 4 lines per inch. If finer ruled paper is used, the 
major division lines should be drawn in with black ink, omit- 
ting the finer divisions. In the case of finely ruled paper, only 
blue-lined paper can be accepted. Tracing paper and blue- 
prints are not acceptable. Lettering and all markings must 
be large enough to be readable after reduction. Mail rolled or 
flat, never fold. Drawings which cannot be reproduced (includ- 
ing pencil drawings) will be returned to the author for redraw- 
ing, thus delaying publication of the paper. Photographs 
should be very distinct and show clear black and white con- 
trasts. They must be on glossy white paper. Avoid round and 
oval photographs. 

CAPTIONS AND LEGENDS: Legends or captions must accom- 
pany each drawing or photograph submitted. If written on the 
drawing or photograph, they should be placed below and well out- 
side the part to be reproduced. It is better to place them on 
separate sheets of paper pasted to the back of the drawings or 
photographs. Each table should have a caption such as Table 
1, Table 2, Table 3, etc. Captions should be complete in them- 
selves so as to make the data intelligible to the reader without 
reference to the text. A duplicate list of captions for figures 
should be included as the last page of the manuscript. Use 
“Fig. 1”? (not Figure 1), Figs. 3 and 4, etc., in both the text 
and the numbering of illustrations. In the text, ‘Eq. (1),” or 
“Eqs. (1) and (2)” are preferable to ‘Equation (1).”” In cap- 
tions and legends, except for ‘‘Fig.’’ and ‘‘Eq.”’ and table head- 
ings, write all words in full; do not abbreviate. Avoid placing 
explanatory written matter in the drawings; it should be in 
the text. 


MATHEMATICAL WorRK: Only the very simplest formulae 
should be typewritten; all others should be very carefully written 
in pen and ink, the writing to be large enough so that ample room 
is provided to mark mathematical matter for the printer. A 
considerable space for marking should be allowed above and be- 
low all equations. All complicated equations should be repeated 
on separate sheets with plenty of space left for marking. The 
solidus should be used for simple fractions appearing within 
the text. Make all expressions clear to the typesetter. Greek 
letters used in formulae should be clearly designated by name 
on the margin of the manuscript. All symbols should be clearly 
written and carefully checked. The difference between capital 
and lower-case letters should be clearly distinguished and care 
taken to avoid confusion between zero (0) and the letter (0), 
between the numeral (one) and the letter (ell) and the prime 
(‘), between alpha and a, kappa and k, u and mu, v and nu, n 
and eta. All subscripts and exponents should be clearly marked 
and dots and bars over letters or mathematical expressions should 
be avoided. Avoid complicated exponents and _ subscripts. 
When it is necessary to repeat a complicated expression, it should 
be represented by some convenient symbol. 


NOMENCLATURE AND ABBREVIATIONS: The National Advisory 
Committee for Aeronautics Nomenclature should be used in pref- 
erence to any others. Standard abbreviations should be used, 
and it should be noted that most abbreviations are lower case, 
such as m.p.h., b.m.e.p., i.hp., b.hp., hp., ... ete. 




















